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ABSTRACT 

Modern large-scale surveys have allowed the identification of large numbers of white 
dwarfs. However, these surveys are subject to complicated target selection algorithms, 
which make it almost impossible to quantify to what extent the observational biases 
affect the observed populations. The LAMOST (Large Sky Area Multi-Object Fiber 
Spectroscopic Telescope) Spectroscopic Survey of the Galactic anti-center (LSS-GAC) 
follows a well-defined set of criteria for selecting targets for observations. This advan¬ 
tage over previous surveys has been fully exploited here to identify a small yet well- 
characterised magnitude-limited sample of hydrogen-rich (DA) white dwarfs. We de¬ 
rive preliminary LSS-GAC DA white dwarf luminosity and mass functions. The space 
density and average formation rate of DA white dwarfs we derive are 0.83±0.16 x 10“^ 
pc“^ and 5.42 ± 0.08 x 10“^^ pc“^ yr~^, respectively. Additionally, using an existing 
Monte Carlo population synthesis code we simulate the population of single DA white 
dwarfs in the Galactic anti-center, under various assumptions. The synthetic popu¬ 
lations are passed through the LSS-GAC selection criteria, taking into account all 
possible observational biases. This allows us to perform a meaningful comparison of 
the observed and simulated distributions. We find that the LSS-GAC set of criteria 
is highly efficient in selecting white dwarfs for spectroscopic observations (80-85 per 
cent) and that, overall, our simulations reproduce well the observed luminosity func¬ 
tion. However, they fail at reproducing an excess of massive white dwarfs present in 
the observed mass function. A plausible explanation for this is that a sizable fraction 
of massive white dwarfs in the Galaxy are the product of white dwarf-white dwarf 
mergers. 

Key words: (stars:) white dwarfs; stars: luminosity function, mass function 


1 INTRODUCTION 

White dwarfs (WD) are the typical endpoint of the evolu¬ 
tion of most main sequence stars. Because nuclear reactions 
do not occur in their deep interiors, the evolution of WDs 
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can be considered as a relativel y simple and well und erstood 
gravothermal cooling process l Althaus et al.ll2010al ). Actu¬ 
ally, the evolutionary co oling times are now very accurate 
(e.g. iRenedo et aLllioid ). providing a reliable way of mea¬ 
suring the WD cooling age from the temperature and sur¬ 
face gravity measured observationally. WDs are hence very 
useful tools in astronomy. For example, the WD luminosity 
function (LF) is an important statistical instrument which 
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has been used not only to derive an accurate age of the 
Galactic disk in the solar neighborhood (W inget e t al.ll 198^ : 
iGarcia-Berro et al.lll98^ : iFontaine et aL T200l|) b ut a lso to 
constrain the local sta r formation ra te ( Noh fc Scalol Il99d : 
IPiaz-Pinto et ah IH; i owe il l2013h . Moreover, the WD 
mass function (MF) has been successfully employed over the 
years as a tool to test the theory of stellar evolution, offering 
information on stellar mass loss. This function has also been 
used to study how the evolution of close b inaries is able to 
produce low-mass WDs (iLiebert et al.llioosi '). thus helping to 
asses the contribution of the distinct evolutionary scenarios 
in producing the current population of Galactic binaries in 
which one of the components is a WD. Finally, the WD age 
function (AF) is also a valuable tool for constraining the 
WD formation history JHu et al. I [2003). 

With the advent of modern, large-s cale surveys such 
as the Sloan Digital Sky Survey (SPSS; Yorl^et^ 2000l l 
or the Super-Gosmos Sky Survey llHamblv et al.ll2001 '). the 
size of the current observational WD samples has increased 
dramatically. This has allowed produ c ing more accurate 
LFs, MF s and AFs ([Ha rris et ajTl2006l : iKepler et al.l 20071 : 
Hu et al.] I2OO7I : iDe Gennaro et ahl 20081 : Rowell fc Hambl-J 


2011^ . However, the major drawback of most studies is the 


complicated target selection algorithms, which incorporate 
observational biases that are almost impossible to quantify 
in numerical simulations that aim at reproducing the ensem¬ 
ble properties of the observed samples. 

Significant observational efforts have also allowed to 
unveil the population of WDs within 20 pc of th e Sun 
llHolberg et al.[l2003 . l2008l : iGiammichele et ^l2012l ~). This 
sample is considerably less numerous than those mentioned 
above, however it is claimed to be reasonable complete and 
can therefore be considered as a volume (rather than mag¬ 
nitude) limited sample, which suffers effectively from no se¬ 
lection biases. 

The LAMOST (Large Sky Area Multi-Object Fiber 
Spectroscopic Telescope) Spectroscopic Survey of the Galac - 
tic anti-center fLSS-GAC: lLiu et al.l]2oi3 : lYuan et al.ll2015l l 
follows a well-defined selection criteria aiming at providing 
spectra for stellar sources of all colours in the Galactic anti- 
center (including WDs) so that they can be studied in a 
statistically meaningful way. LSS-GAG started operations in 
2011 and will provide a signihcantly larger sample of WDs 
than those within 20 pc of the Sun. In this paper we derive 
preliminary observed LF, MF and AF of WDs identified 
within the first data release of LSS-GAC, and use a state- 
of-the-art Monte Carlo population synthesis code adapted to 
the characteristics of the survey to simulate the WD popu¬ 
lation in the Galactic anti-center. We apply the LSS-GAC 
selection criteria to the simulated samples, carefully eval¬ 
uate all possible observational biases, and derive synthetic 
LFs, MFs and AFs. This exercise allows us to perform a 
meaningful comparison between the outcome of simulations 
and the observational data. 


2 THE LAMOST SPECTROSCOPIC SURVEY 
OF THE GALACTIC ANTI-CENTER 

LAMOST is a quasi-meridian reflecting Schmidt telescope 
located at Xinglong Observ i ng Station in the Hebei province 
of China (ICui et al.1 12OI2I : ILuo et ^ I2OI2I I . The effective 


aperture of LAMOST is about 4 meters. LAMOST is ex¬ 
clusively dedicated to obtain optical spectroscopy of celes¬ 
tial objects. Each “spectral plate” refers to a focal surface 
with 4,000 auto-positioned optical fibers to observe spectro¬ 
scopic plus calibration targets simultaneously, equally dis¬ 
tributed among 16 hber-fed spectrographs. Each spectro¬ 
graph is equipped with two GCD cameras of blue and red 
channels that simultaneously provide blue and red spectra 
of the 4,000 selected targets, respectively. 

The LSS-G AC is a majo r cor nponent of t he LAMOST 
Galactic survey (iDeng et al.ll20l3 : IZhao et al.l[20l3l . By se¬ 
lecting targets uniformly and randomly in (r, g — r) and 
(r, r — i) Hess diagrams llYuan et al.ll2015ll . the LSS-GAG 
main survey aims at collecting AA3,700 - 9,000 low resolu¬ 
tion {R ~ 1, 800) spectra for a statistically complete sample 
of ~ 3 million stars of all colours down to a limiting magni¬ 
tude of r =17.8 mag (18.5 mag for limited fields), distributed 
in a contiguous sky area of over 3,400 deg^ centered on the 
Galactic anti-center (|&| ^ 30°, 150 ^ I ^ 210°). The simple 
yet non-trivial target selection of the survey makes possi¬ 
ble to study the underlying stellar populations for any given 
type of target, such as WDs. In addition to the main survey, 
the LSS-GAG also includes a survey of the M 31/M 33 area, 
targeting hundreds of thousands of objects in the vicinity 
helds of M31 and M33, and a survey of Very Bright (VB) 
plates, targeting over a million of randomly selected very 
bright stars (r ^ 14 mag) in the northern hemisphere dur¬ 
ing bright/grey nights. 

Targets for the LSS-GAC main survey are selected for 
spectroscopic follow-up from the Xuyi Schmidt Telescope 
Photometric Sur vey of the Galactic anti-center (XSTPS - 
GAG) catalogue llLiu et al.l l20l4 IZhang et al.1 l2013l . l2014l 'l . 
The XSTPS-GAG survey was initiated in October 2009 and 
completed in March 2011. Photometry in SDSS g, r and 
i bands was acquired with the Xuyi 1.04/1.20 m Schmidt 
Telescope located at the Xuyi station of the Purple Moun¬ 
tain Observatory, P. R. China. The XSTPS-GAC surveyed 
an area of about 5,400 deg^, from RA ~ 3 to O** and 
Dec ~ —10 to -1-60° to cover the LSS-GAC main survey 
footprint, plus an extension of ~ 900 deg^ to the M31, M33 
region to cover the LSS-GAC M 31/M 33 survey area. The 
resulting catalogues contain about 100 million stars down to 
a limiting magnitude of r ~ 19.0 mag (l Ocr), with an astro - 
metric calibration accuracy of 0.1 arcsec (IZhang et al.ll20i3l 
and a global pho tometric calibration accuracy of ~2 per cent 
llLiu et al.ll2014l . Yuan et ah, in prep.). 

For the LSS-GAC main survey, bright (B), medium- 
bright (M) and faint (F) plates are designated to tar¬ 
get sources of brightness 14.0< r < 16.3 mag, 16.3 < 
r <17.8 mag and 17.8 < r ^ 18.5 mag, respectively. During 
the Pilot (Oct. 2011 - Jun. 2012) and the hrst year Regu¬ 
lar (Oct. 2012 - Jun. 2013) Surveys of LAMOST, a total of 
1,042,586 [750,867] spectra of a signal to noise ratio (S/N) 
10 at 7450 A [S/N ^ 10 at 4650A] have been collected, 
including 439,560 [225,522] spectra from the main survey. 
Most of the stars are from the B plates. 

The raw data were reduced with the LAMOST 2D 
pipeline (Version 2.6) fe.g. ILuo et al.|[2004h . following the 
standard procedures of bias subtraction, cosmic-ray removal, 
ID spectral extraction, flat-fielding, wavelength calibration, 
and sky subtraction. The LAMOST spectra are recorded in 
two arms, 3,700 - 5,900 A in the blue and 5,700 - 9,000 A in 


© 2015 RAS, MNRAS 000, ??-?? 

































































DA white dwarfs from the LSS-GAC 3 







Figure 1. Reduced Xspec ~ S/N distribution obtained by fitting 
one of our WD templates [Tgff=14,000K, log(g)=7 dex, see Fig.[2] 
to the entire LSS-GAC spectroscopic data base. Objects falling 
in the area defined by Xspec < 0.002 X S/N^'^ (red curve) are 
considered WD candidates (magenta dots). 

the red. The blue- and red-arm spectra are processed in¬ 
dependently in the 2D pipeline and joined together after 
flux calibration. Flux calibration for the LSS-GAC observa¬ 
tions is performed usi ng an iterative algo rithm developed at 
Peking University by IXiang et al.l (l2015l ). achieving an ac¬ 
curacy of about 10 per cent for the whole wavelength ranges 
of blue- and red-arm spectra. A detailed description of the 
target selection algorithm, survey design, observations, data 
reduction, and value-added catalogs of the hr st LSS-GAC 
data release is presented in lYuan et al.l ll2015h . In the cur¬ 
rent work we only consider spectra obtained by the main 
LSS-GAC survey. A list including the names of the plates 
employed for the main survey is provided in Table A1 (see 
Appendix). 


3 THE LSS-GAC WD SAMPLE 

In this section we describe our methodology for identifying 
WDs within the LSS-GAC spectroscopic data base. We do 
this following two independent but complementary routines. 
The first identifies WDs by x^-teniplate fitting all LSS-GAC 
spectra, the second selects WDs by applying a well-defined 
colour cut to the XSTPS-GAC photometric catalogue. We 
also estimate the spectroscopic completeness of the LSS- 
GAC WD sample. 

3.1 The fitting method 

We use the x^-template f itting method described 
by iRebassa-Mansergas et al.l (I201C ) to identify LSS- 
GAC WDs. This routine was originally developed to 
search for SPSS WD-main s e quence bina ries (WDMS; 
IRebassa-Mansergas et aki l201Cll . l2012l . l2013l) . However its 
mathematical prescription can be easily implemented to the 
LSS-GAC data for identifying single WDs. As the first step, 
a given WD template is used to fit the entire LSS-GAC 
spectroscopic data base and the resultant reduced Xspec 
values are recorded. Those Xspec > together with the overall 


Figure 2. Four WD template examples used in this work. Ef¬ 
fective temperature and surface gravity are indicated for each of 
them. 

S/N ratios of the LSS-GAC spectra are represented in a 2D 
map and an equation of the form, 

Xma. = a X SN^ (1) 

is defined so that all systems below this curve (Xspec < Xmax) 
are considered as WD candidates, where a and b are free 
parameters dehned for each template. An example is illus¬ 
trated in Fig.[T] The form of Eq. (1) is dehned to account 
for systematic errors between template and observed spec¬ 
tra, which become more important the larger the S/N is. 
The particular shape of the curve is carefully evaluated for 
each template as a compromise between the number of ex¬ 
cluded/selected spectra. In other words, we aim at selecting 
as many WD candidate spectra as possible of similar spec¬ 
tral features as the templates used (note that, strictly speak¬ 
ing, the only free parameter we are considering in the htting 
process is the normalization factor). However, the number 
of selected targets cannot be too large, as otherwise the au¬ 
tomatization becomes useless. 

The spectra of all WD selected candidates are then vi¬ 
sually inspected and systems that are not WDs are simply 
excluded. This exercise is repeated for all WD templates con¬ 
sidered. In this case, our templates are the result of adding 
artihcial Gaussian noise to 45 DA (hyd rogen-rich) car efully 
selected WD model spectra provided bv lKoest^ (l2010l ) that 
cover a broad range of WD effective temperatures (6,000K- 
60,000 K) and surface gravities (7 < log g <9 dex). The 
model spectra were binned to the LSS-GAG resolving power. 
Four of our considered WD templates are shown in Fig. [21 
Obviously, using only DA WD templates restricts our search 
to hydrogen-rich WDs. However, as we show below (Sec- 
tion l3.3ll . the number of non-DA WDs currently observed 
by the LSS-GAC is rather small. Moreover, it has to be 
stressed that the presence of a small fraction of non-DA 
WDs d oes not seem to influe nce significantly the shape of 
the LF llGoiocaru et al.l[2014h . 

As already mentioned, LSS-GAC spectra are the result 
of combining a blue and a red optical spectrum obtained 
from two separate arms. Given that single WDs are blue ob¬ 
jects, in principle, the blue arm spectra alone should allow 
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r 1 mag 

Figure 3. (g — r) versus (r — i) colour-colour plot illustrating 
the location of LSS-GAC DA white dwarf candidates (black solid 
dots within the colour cut defined by the red solid lines). The 
92 confirmed LSS-GAC DAs are shown as cyan dots, others as 
grey dots The g magnitudes of 22 LSS-GAC DA white dwarfs 
are found to be unreliable, which explains why these WDs are 
found outside the expected colour region (two of them beyond 
the displayed area). 


US to efficiently identify snch stars. To investigate this, we 
separately applied the above described rontine to all blue- 
arm LSS-GAC spectra as well as to all combined (blue- and 
red-arm) spectra and compared the results. We found that 
the number of identified WDs differed significantly. Specif¬ 
ically, we identified ~ 30 per cent more objects using the 
combined, blue plus red, spectra. Moreover, all WDs iden¬ 
tified using only the blue-arm spectra were included in the 
list of targets found using the combined LSS-GAC spectra. 
The overall shape of the continuum spectrum is hence re¬ 
quired for efficiently selecting single WDs. We found this ef¬ 
fect to be most important when selecting WDs against hot 
(A) main sequence stars in which the strengths of the Balmer 
lines are similar. This was also the case when selecting WDs 
amongst low S/N (< 5) spectra. For these two reasons, we 
decided to only consider the combined LSS-GAC spectra 
of S/N ratio above 5 in this work. This excluded 15,442 
LSS-GAC spectra for which the default 2D pipeline (Sec- 
tion[2]) failed to deliver usable red- and/or blue-arm spectra. 
It also excluded 622,112 spectra of S/N ratios below 5. This 
requirement is further supported by the fact that reliable 
stellar parameters (necessary for obtaining reliable LF, MF 
and AF) cannot be derived for low S/N ratio WD spectra 
llRebassa-Mansergas et al.ll201Cll ). Thus, the total number of 
LSS-GAC spectra considered in this work was 306,600, of 
which 20,029 were selected as WD candidates by our x^~ 
template fitting method. A visual inspection revealed that 
94 of those are genuine WDs, 81 of DA type (see Table[T]). 
We also identified 26 WDMS binaries. With the exception of 
J080357.21-1-334138.6, these binaries are al ready included in 
the LAMOST WDMS binary catalogue by I Ren et al.l ll2013l . 
l2014h . 


Table 1. Number of WDs (including DAs, DBs, etc.) and DAs 
alone identified by the y^-S/N template fitting and colour meth¬ 
ods, in common by the two methods, independently by each 
method, and missed by both methods. The total number of LSS- 
GAG WDs and DAs are provided in the last row. 



Awd 

Nda 

X^-S/N method 

94 

81 

Colour method 

78 

70 

Common by both 

63 

57 

X^-S/N method alone 

27 

20 

Colour method alone 

15 

13 

Missed by both 

2 

2 

Total 

107 

92 


3.2 The colour cut method 

In the previous section we searched for DA WDs by x^~ 
template fitting all LSS-GAC spectra. Here we complement 
this with an independent strategy based on XSTPS-GAC g, 
r and i photometry (Section[2} . This exercise relies simply on 
applying a well-defined cut to the XSTPS- GAC g—r and r—i 
colours (see Tablet o fj Girven et al.ll20li] and Fig.O of this 
paper). iGirven et al.l (1201 1[) provide two additional colour 
cuts, using the u and z magnitudes respectively, however 
these are not used here because the XSTPS-GAC does not 
provide magnitudes in those filters. 16,824 sources fall within 
the area defined by our cut and visual inspection of spectra 
of those targets confirms 78 as WDs, of which 70 are DAs 
— see Table[T] The vast majority of the remaining targets 
are single A, F main sequence stars whose colours overlap 
with those of cool WDs. The contamination by quasars is 
negligible as these sources are generaly too faint to be ob¬ 
served by the LSS-GAC survey. 63 (57) of these 78 (70) WDs 
(DAs) are already discovered following the above described 
X^-template fitting method. Thus, the colour method adds 
15 new objects (13 of which are DAs) to our sample. The to¬ 
tal number of LSS-GAC WDs thus raises to 105, and among 
them we classify 90 as DA WDs (Table[T]). The complete 
LSS-GAC DA WD sample is provided in Table A2 (see Ap¬ 
pendix) . 

Two important conclusions can be drawn from the 
above exercise. First, the x^-template fitting method failed 
to identify ~ 15 per cent of the whole LSS-GAC DA WD 
sample. Visual inspection of the spectra of those systems 
revealed that they were either of low S/N ratios (~ 5 — 6 
in the blue-arm spectra), or were subject to artifacts of a 
bad flux calibration/merging of the two-arm (blue plus red) 
spectra, or a combination of both. Although those effects 
clearly affect the identification of DA WDs when applying 
the x^-template fitting method, they do not influence the 
colour selection of DA WDs. Secondly, the DA WD colour 
cut we applied missed ~ 22 per cent of the whole DA WD 
sample identified. A closer inspection revealed the g mag¬ 
nitudes of these targets to be unreliable. This caused those 
objects to fall far from the colour locus expected for DA 
WDs (see Fig. EJ. 
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Figure 4. Spectral model fits to two DA WDs in our sample. Top row, 1st and 3rd panels: best-fit (black lines) to the normalised H/3 
to He (grey lines, top to bottom) line profiles. Top row, 2nd and J^th panels: 1, 3 and 6it contour plots in the Tgfj — logg plane. The 
black contours refer to the best line profile fit, the red contours to the fit of the entire (continuum plus lines) spectrum. The dashed 
line indicates the occurrence of maximum H/3 equivalent width. The best “hot” and “cold” line profile solutions are indicated by black 
dots, the best fit to the entire spectrum (continuum plus lines) is indicated by a red dot. Bottom row, upper panels: the WD spectra and 
associated uncertaintied (grey lines) along with the best-fit model spectrum (black line) for the 4,000—5,500 A wavelength range. The 
that result from the Balmer line fitting (cold and hot solutions) and the entire spectrum fitting are also indicated. Bottom row, lower 
panels. The residuals of the fit in the 4,000-5,500 A wavelength range. The Tgff and logg values listed in Table A2 are determined from 
the best Balmer line profile fit. The fit to the entire spectrum is only used to select between the “hot” and “cold” solutions of Balmer 
line profile fit (see also Fig.lSt. 


3.3 Spectroscopic completeness of the LSS-GAC 
DA WD sample 

The spectroscopic completeness is defined as the fraction of 
DA WDs that we have identified compared to the total num¬ 
ber of WDs observed by the LSS-GAC. It can be assumed 
that the DA WD spectroscopic sample is 100 per cent com¬ 
plete within the colour selection box defined by the g — r 
and r — i colours, as we have visually inspected every sin¬ 
gle spectrum within this region. There are 70 DA WDs in 
this colour box, of which 57 were found by the y^-template 
fitting method (Table[T]). The spectroscopic completeness of 
the method within the colour box is therefore 80 per cent. 
The number of DA WDs found by the ^(/^-template fitting 
method outside the colour box is 20 (Table[T|. If we assume 
that the spectroscopic completeness of the x^-template fit¬ 
ting method is not strongly colour dependent, then the total 
number of DA WDs we expect to lie outside the colour box 
is 20/0.8 = 25. The total number of DA WDs that could 
have been detected is therefore (70-1-25) = 95, and the over¬ 
all spectroscopic completeness is (70-l-20)/95 ~ 95 per cent. 
The two independent methods we have used to identify LSS- 
GAC DA WDs, namely the x^-template fitting method and 
the colour cut method, seem to complement each other very 
well and manage to identify the vast majority of DA WDs 
observed by the LSS-GAC survey. 

An addional way to quantify the spectroscopic com¬ 
pleteness is by cross-correlating our list with the LAM- 


OST DA WD ca talogu es published by IZhao et~n! (201^ 
and IZhang et al.l (|201 ^. Of th e 16 DA WDs of I Zhao et ai.l 
ll2013l) . and 27 of Zhang et al.l (l2013l) that fall in the area 
observed by LSS-GAC and have spectra of S/N ratios ^ 5 
in both the red- and the blue-arm spectra, 15 and 26 are 
within our sample, respectivelj0. Therefore, we have missed 
two DA WDs (Table[T]). This exercise suggests a spectro¬ 
scopic completeness of 94 — 96 per cent for our sample, in 
good agreement with the above estimated value. 

In Table A2 we include the two DA WDs that we have 
missed. This increases our LSS-GAC WD catalogue to 107 
WDs. Of these, 92 are of the DA type (Table[T]). Given that 
DA WDs ar e the most common o nes - ~85 per cent, see 
for example iKleinman et al.l (l2013l ) - and our final sample 
contains 92 DA WDs and seems to be highly complete, it can 
be said that the number of non-DA WDs currently observed 
by LSS-GAC is very small. 


4 STELLAR PARAMETERS AND DISTANCES 

We determine the s tellar parameters following the fitting 
routine described bv iRebassa-Mansergas et al.l ll2007l ). This 


^ In passing we note that the total number of objects from the 
catalogue bv IZhang et al.l ll2013l l that fall within the Galactic anti- 
center area is 71, however we exclude 44 as they are not WDs. 
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Figure 5. . Observed LSS-GAC spectra of the two DA WDs shown in Fig.[4| (gray), together with the best Balmer line model fits 
(black solid line for the hot solutions, black dashed line for the cold solutions). The GALEX near- and far-ultraviolet fluxes and the 
XSTPS-GAC optical fluxes (derived from the p,r, 2 magnitudes) are shown as red solid dots. Left panel: the predicted ultraviolet and 
optical fluxes clearly select the hot solution of the Balmer fitting, in agreement with the solution preferred by fitting the entire LSS-GAC 
spectrum (Fig. 01 top row, second panel). Right panel: the GALEX ultraviolet fluxes are in good agreement with the cold Balmer line fit 
solution and thus clearly rule out the hot solution favoured by the fit to the entire LSS-GAC spectrum (Fig.[4] top row fourth panel). 
Note that in this particular case the g XSTPS-GAC magnitude is found to be unreliable (Fig. 0)1. 


relies on using the entire model grid of DA WDs of lKoest^ 
1 I 2010 I) to fit the normalised H/3 to He line profiles of each 
WD spectrum for determining the effective temperature 
(Tefi) and surface gravity (logg). Two examples are illus¬ 
trated in Fig. |4] (top row, first and third panels). Given that 
the equivalent widths of the Balmer lines go through a max¬ 
imum near T^s = 13, 000 K (with the exact value being a 
function of logg), Teff and logg determined by fitting the 
Balmer line profiles are subject to an ambiguity, often re¬ 
ferred to as “hot” and “cold” solutions. This implies that, 
at a given equivalent width, the Balmer lines of a cold and 
massive WD have the same profile as the Balmer lines of 
a hot and less massive WD (see Fig.jd] top row second and 
fourth panels). We break this degeneracy as follows. First, 
we fit the entire WD spectrum (continuum plus lines) with 
the same grid of model spectra used for fitting the Balmer 
lines (see bottom panels of Fig. |4]). We used the 4,000-5,500 A 
wavelength range covered by the blue arm of the LSS-GAC 
spectra. This is because in some cases the combined LSS- 
GAC spectra are subject to artifacts of a bad merging of the 
blue plus red arm individual spectra. Since the continuum 
spectrum of a DA WD is mostly sensitive to Teg, the best-fit 
value from the entire spectrum indicates which of the two 
solutions is the prefered one. However, because of uncertain¬ 
ties in the flux calibration, the fit to the entire spectrum can 
only be used for breaking the degeneracy between the hot 
and cold solutions, rather than for obtaining a reliable set 
of stellar parameters. This in turn implies that the solution 
prefered by the best-fit to the entire spectrum may be sub¬ 
ject to systemmatic uncertainties. Thus, in a second step 
the choice between hot and cold solution is further guided 
by co mparing the ultraviolet GALEX (Galaxy E volution Ex¬ 
plorer; |MartiFeF^i][20^;[Mprrisse^^2laI][200^ and optical 


XSTPS-GAC0 fluxes to the fluxes predicted from each so¬ 
lution (see two examples in Fig.[5]). This exercise allows us 
to confidently select the correct solution for 75 of our 92 
LSS-GAC DA white dwarfs. We use these 75 LSS-GAC DA 
white dwarfs as the sample of analysis in this work. 

ft has been shown that spectroscopic fits that use ID 
model atmosphere spectra such as those employed in this 
work result in systematically overesti mated surface gravi- 
ties for WDs cooler than ~12,000K (iKoester et al.l 120091 : 
iTremblav et al.l 1201 ill. We have t hus applied the 3D cor¬ 
rections of Tremblay et al.l (l2013li to Teg and loggr deter¬ 
mined above. We t hen interpolated the Teg and logg val¬ 
ues in t he tables ofiRen edo et al.l ll2010 li and AiHmus et al.l 
ll2005ll : I Althaus et al.l 1 2007 1: Althaus et al.l (2010 lJ) to ob¬ 
tain masses, cooling ages, absolute Mg, Mr and Mi and 
bolometric (Mboi) magnitudes for our WDs. These cooling 
sequences provide absolute magnitudes in the UBVRI sys¬ 
tem, which are converted into the ugriz system using the 
equations of Ijordi et al.l ll2006li . Distances were finally ob¬ 
tained from the distance moduli of our targets, with the 
XSTPS-GAC (/,r, i magnitudes corrected for ex tinction us¬ 
ing th e 3D Galactic extinction map provided bv IChen et al.l 
ll2014|) . 


The effective temperature, surface gravity, mass, cool¬ 
ing age and distance distributions of the 75 DA WDs for 
which we are able to break the degeneracy between the hot 
and cold solutions obtained by fitting the Balmer lines of 
their LSS-GAC spectra are shown in Fig.|^ Inspection of 
this figure reveals that the vast majority of LSS-GAC WDs 
are located at short distances (~50-300pc). The mass and 
log g distributions broadly peak at 0.6 Mq and log p = 8 dex 


^ The optical fluxes are derived directly from the XSTPS-GAC 
g, r, i magnitudes. In those cases in which the g magnitudes are 
unreliable (Fig.[3]l we either use only the r, i fluxes as guidance 
or substitute the g magnitudes by those provided by SDSS, when 
available. 
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Figure 7. Correlations between the effective temperatures, masses, distances and cooling ages of the 75 LSS-GAC DA white dwarfs 
shown in Fig. 0 (black dots). We also show a typical realization of the Monte-Carlo simulations described in SectionO(red stars). As can 
be seen the degree of overlap of the synthetic sample and the observational one is nearly perfect in all six panels. 


respectively, features th at have been observed previously in 
numerous studies ( e.g. [ Koester et al.l 1 19791 : iHolberg et al.l 
I2OO8I : iKepler et al.l l201a) . The effective temperatures are 
clustered between ~ 10,000-15,000 K, and the majority of 
WDs have cooling ages between 0.01 and 1 Gyr. In Fig.[7]we 
show a grid of panels displaying the correlations between the 
effective temperatures, masses, distances and cooling ages 
(black dots) . We do not include the surface gravity here, 
as it is nearly equivalent to mass. It becomes clear that the 
higher the WD effective temperature the larger trends to be 
the distance, a clear selection effect indicating that cooler 
WDs need to be generally closer to us to be observed by 
the LSS-GAC. Also, hotter WDs trend to be further away, 
as otherwise they saturate the lower magnitude limit of the 
LSS-GAC. These effects are also noticeable when inspect¬ 
ing the relation between mass and distance: low-mass WDs 
are brighter and trend to be found at larger distances. It 
also becomes evident that, as a simple consequence of the 
WD cooling, the cooling ages decrease for increasing effective 
temperatures. Because of this there is also a tight relation 
between cooling age and distance, i.e. shorter cooling ages 
imply hotter WDs, therefore further distances. There is also 


a clear correlation between the cooling age and the distance, 
a quite natural behavior, since given that the mass distribu¬ 
tion of WDs has a narrow peak around 0.6 Mq, the cooling 
ages and the effective temperature are very tightly correlated 
(see the leftmost central panel of this figure). Also shown in 
this figure is a typical realization of the Monte Garlo simu¬ 
lations described below (Section [G]). It is quite apparent the 
high degree of overlap between the synthetic and the ob¬ 
served samples, indicating that our models reproduce with 
high fidelity the selection procedures, and the astronomical 
properties of the observed population. 


From the above discussion it becomes clear that our 
observational sample and corresponding distributions are af¬ 
fected by selection effects typical of those incorporated by 
magnitude limited surveys. We correct for those observa¬ 
tional biases in the next section. 
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the LSS-GAC survey. This requires considering the lower 
and upper magnitude limits of each of the 16 spectrographs 
of each LSS-GAC plate. For each spectropgraph, the lower 
and upper magnitude limits dehne respectively the mini¬ 
mum (dmin) and maximum (dmax) distances (and therefore 
minimum and maximum volumes, Vmin and Vnia.A at which 
the considered WD would have been detected. The total 
maximum volume of a WD, Vwd, is the sum over the indi¬ 
vidual m aximum volumes obtained from each spectrog raph 
- see also IHu et al.l l^2007^ ; |Limoges fc Berger^ ll201ofl : 


^spec 

^ . r“max 

OJi ' 


Vwd = Vmax - Vmin = ^ ^ ^ AiYv'^dr = 
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Figure 6. From top to bottom: mass, effective temperature, sur¬ 
face gravity, cooling age and distance distributions of the 75 LSS- 
GAC DA white dwarfs for which we are able to break the degen¬ 
eracy between the hot and cold solutions obtained from fitting 
the Balmer lines of their spectra. 


5 THE LSS-GAC DA WD LUMINOSITY AND 
MASS FUNCTIONS, AND THE DA WD 
FORMATION RATE 

Magnitude limited surveys such as the LSS-GAC survey are 
affected by selection effects. Therefore, any parameter dis¬ 
tribution that results from the analysis of a given observed 
population is subject to observation al bia ses. The 1/Um ax 
method described in ISchmidj lll968h and iGreenI lIlQSd j is 
aimed at removing these biases. In our case this is done cal¬ 
culating the maximum volume in which each of our WDs 
would have been detected given the magnitude limits of 


where b is the Galactic latitude of the WD, and toi is the 
solid angle in steradians covered by each spectrograph (1.2 
deg^ X TT^/ISO^; 1® ^li® total area observed by the 

survey, also in steradians). The factor takes into ac¬ 

count the non-uniform distribution of stars in th e direc¬ 
tion perpendicular to the Galactic disc llFeltenlll976ll . where 
z = r X sin(6) is the distance of the WD from the Galac¬ 
tic plane, and zp is the scale height, which is assumed to 
be 250 pc llLiebert et ahll^OSl : IHu et aDl2007ll . In the cases 
where two or more spectrographs observe the same region of 
sky, we consider the overlapping region with the largest vol¬ 
ume, spanning between the smallest lower magnitude limit 
and the highest upper magnitude limit of the overlapping 
spectrographs. 

Once Vwd is calculated for each WD in the observed 
sample the space density of WDs is simply obtained as 
X) 1/VwD = 0.83 ±0.16 X 10“® pc~®, where the sumation 
is over all the WDs in the sampl^fj. However, it has to be 
noted that the space density derived here represents an ab¬ 
solute lower limit, as we are able to derive reliable stellar 
parameters for only 75 of the 92 DAs in our sample. That 
is, we are considering just 81 per cent of the observed sample 
in the analysis. Moreover, the lowest effective temperature 
value among LSS-GAG DA white dwarfs is ~6,500 K. WDs 
of lower effective temperatures are too faint to be detected 
by the survey, and are therefore not accounted for in our cal¬ 
culation of the space density. The space density as a function 
of the bolometric magnitude Mboi, mass Mwd and cooling 
age tc yield the WD LF, MF and AF, respectively. Each of 
these functions is analysed in the following sub-sections. 

The 1 /Umax method described above can be also used to 
quantify the completeness of the observed sample, i.e. the 
percentage of WDs that are still missing because of selec¬ 
tion effects after applying the 1/Vmax method. This com¬ 
pleteness must not be taken as the spectroscopic complete¬ 
ness of the LSS-GAG sample, which is 95 per cent (see Sec- 
tion l3.3l) . If the sample is complete, th en the avera ge value 
{V - Vmin)/(Umax “ Vmin) should be 0.5 (iGreeijl 1980l) (where 


^ The error of the space density was obtained artificially produc¬ 
ing 200 versions of the observed luminosity function by varying 
the value of the function in each bin with a random value sam¬ 
pled from a poisson distribution proportional to the error bar 
corresponding to that bin. 
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Figure 8. The LF of LSS-GAC DA white dwarfs and associated 
errors (solid black line and solid do ts). For the sake of compa rison 
we also show the LFs obtai ned bviDe Genna ro et al.l ll 20081) (red 
dashed line and red st ars). [ T orres et^L ll20 14) (blue dotted line 
and blue squares), and lGiarmnichele et al.l i|2012l) (solid grey line 
and grey triangles). 


V is the volume of the WD as defined by its distance, i.e. the 
same as Equation but integrating from 0 to d). In our 
case this quantity is 0.4, which corresponds to a complete¬ 
ness of 80 per cent. Of course, the above derived estimate 
of the completeness is within the context of the magnitude 
limits of the LSS-GAC survey, i.e. it does not account for 
populations of WDs that are too faint/rare to make it into 
the observed sample. Moreover, 19 per cent of the observed 
sample has not been considered in the analysis. If we were 
able to constrain the stellar parameters of these WDs, then 
the completeness would increase. 


Figure 9. The MF of LSS-GAC DA white dwarfs. 


Inspection of Fig. [8] reveals that, for Mboi 6 mag, the 
L F derived in this work is in good agreement with the LF 
oflpT Gennaro et al.l (|2nn8ll . The appare nt disagreement be¬ 
tween our LF and the one obtained bv iGiammichele et al.l 
ll2012ll is likely due to the fact that the latter study includes 
all WDs (not only DAs) in a volume-limited local sample 
which does not presumably suffer from completeness issues, 
and therefore the space density is higher. For Mboi < 6 mag, 
our LF is also in broad agreement with that of iTorres et al.l 
ll2014ll for hot DAs. However, the number of LSS-GAC WDs 
falling in these bins is too small for a meaningful compar¬ 
ison between the two studies. It should be also noted that 
the WD LF actually continues to fainter magnitudes than 
those shown in Fig.O however we do not display those bins 
as these objects are too faint to be present in the LSS-GAC 
sample. 


5.1 The luminosity function 

The LSS-GAC WD LF is shown in Fig.El its associated er¬ 
rors are calculated following iBovld (ll989lL For compar ison 
we show the LFs obtained by De Gennaro et al.l (l2008li for 
the SDSS survey, [Torres et al.l 1 20141') for hot DAs in the 
SPS S (which supersede s the LF of Krzesinski et al.ll2009ll . 
and IGiammichele et al.l (l2012l ~l for a local and volume lim¬ 
ited sample of WDs. We note that two more LF s are a vail- 
able from the SDSS, provided by [Harris et al. | (l200d) for 
photometric WD candidates, and by [Hu et all ([20071 1 for 


Data Release 1 spectra of DA WDs. Howeve r, we do not 
show them in Fig. [8] because (1) the LF of [Harris et al.[ 
|[2006ll assumes a constant log g = 8dex and very much 
resembles the one provided by iDe Gennaro et al.l ([20081 1 
when co nsidering DA WDs of th e same logg value; (2) the 
work by pe_Gennaro et ^ ([2008f l supersedes the analysis of 
[Hu et an ~(2007 ). Moreover, we decide not to show the LFs 
obtaine d by Liebert et all ([200511 for the Palomar Green Sur¬ 
vey and lLimoges fc Bergeron (l2010|l for the KISO surve y, as 
they also resemble the LF of iDe Gennaro et al.l ([2008fl but 
contain considerably fewe r objects. Because o f com pleteness 
issues, the LF derived bv [Rowell fc Hamblvl ((^oTJ) for WD 
candidates in the Super-Cosmos survey is not included nei¬ 
ther. To avoid clustering of data in Fig.^we have also opted 
not to show the errors of all the above mentioned LFs. 


5.2 The mass function 

The MF of the LSS-GAC WDs is displayed in Fig.O As ex¬ 
pected, it shows a clear peak around 0.6 Mq. A relatively 
high percentage (~10 per cent) of low-mass (<0.5 M©) WDs 
can also be seen. Traditionally, the existence of these low- 
mass WDs has bee n attributed to binary interactions (e.g. 
[Liebert et ahll^OSll . and indeed it has been demonstrated 
that the majority of low-mass WDs are formed in binaries 
dRebassa-Mansergas et ahl [201 ll : [Kilic etTI [2012[ l. There¬ 
fore, DA WDs in those low-mass bins are expected to be 
part of binaries that contain unseen companions. Inspection 
of Fig. [^ also reveals a large fracion (~30 per cent) of mas¬ 
sive (^0.8 Mq) WDs. a high-mass feature has also been 
regularly detected in a number of studies fe.g. IPiebert et al.[ 
[2(3(351 : [Kepler et al.[ [20071 : [Kleinman et ahlkoisli and it has 
been claimed that it arises as a consequence of WD binary 
mergers. Population synthesis studies however do not pre¬ 
dict more than ~10 per cent of t he entire WD population be¬ 
ing the resu l t of b i nary mergers ( Han^et_^}ll994 JHmIIiOO^ 
Meng et al.[ [20081 : [Toonen et al 20121 : Garcia-Berro et all 


mg 

20121) . Alternatively, a large number of high-mass WDs in 


the MF presented here may be the consequence of large un¬ 
certainties in the mass determinations (the mass errors of 
some of our objects are estimated to be > 0.1 — 0.15M©, 
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which can move some objects across the corresponding mass 
bins). We will further discuss the large percentage of high- 
mass WDs identified in our MF in Section lTTll 


5.3 The formation rate of DA WDs 


We now estimate the average DA WD f ormat ion rate fol¬ 
lowing the method outlined bv iHu et ahl ll2007ll . If the WD 
formation rate is assumed to be constant during the last 
Gyr, then the slope of the cumulative AF (see Fig. HOD can 
be considered as the average formation rate. Thus, we simply 
fit the cumulative AF with a straight line (the red dashed 
line in Fig. fTUll . and identify the slope of 5.42 ± 0.08 x 10“'^® 
pc“® yr~^ of the fit as the DA WD formation rate. Inspec¬ 
tion of Fig. [in] reveals that, as expected, the maximum of the 
cumulative AF is 0.83 x 10~^pc~®, i.e. the total LSS-GAC 
DA WD space density. 

Numerous studies in the past two decades have ob- 
tained WD formati on rates. The most recent analysis 
dVerbeek et al.l [ioi^ results in a birth rate of 5.4 ± 1.5 x 
10 ~^® pc~® in excellent agreement with the value 

estima ted here. The formation rates derived by I Hu et al.1 
ll2007l. 2.5 - 2.7 x 10 "^^ pc"^ yr"^) iLiebert et al.1 (l20^ 
6 X 10“^^ pc“® vr"^!. lilolberg et ^ (l2002l . 6 x 10~^^ pc“® 
yr“^) are also broadly consistent with the value estimated 
here. For comparison, earlier studies yield for mation rates 
that are generally consi derably higher (e .g. iGree iil Il980l. 
20±1 0xl0~” pc~^ yr~MWeidemannlll99ll. 23x10”^® pc"® 
vr~^: IWnnes et aTlll997l . 8.5±1.5 x 10"^® pc"® yr"^). This 
may be a consequence of the recent improvement in qual¬ 
ity and size of WD data sets. Planeta ry nebulae birth rates 
are a lso found to be higher i n genera l Ishida &: Weinberger! 
119871. 80 X 10"^® pc"® vr"^: |PhillipEll200l 21 x 10"” pc~® 
vr"^: [hY^l2008l . 8 ± 3 xlO"^® pc"® yr"®1. T he discrepancy 
has been discussed in detail in lLiebert et al.l (120051 1. 


6 THE LSS-GAC SIMULATED DA WD 

LUMINOSITY, MASS AND CUMULATIVE 
AGE FUNCTIONS 

In the previous sections we presented and characterised the 
sample of DA WDs identified within the data release 1 of 
the LSS-GAC. We also derived the DA WD space density, 
which has been used to construct the preliminary LF and 
MF of LSS-GAC DA WDs. Finally, we estimated the average 
DA WD formation rate from the DA WD cumulative AF. 
In this Section we simulate the LSS-GAC DA WD popula¬ 
tion and take advantage of the well-defined selection criteria 
employed by the LSS-GAC survey to evaluate the fraction 
of simulated WDs that would have been observed by the 
LSS-GAC survey. This will allow us to directly compare the 
ensemble properties of the observational data sets with the 
outcome of the simulations (see Section[7l). 


6.1 The population synthesis code 

We provide here a brief description of our Monte Carlo WD 
population synthesis code, and of the changes done to adapt 
it to reproduce the LSS-GAC. A more complete and de¬ 
tailed description of the principal components of this method 
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Figure 10. The LSS-GAC DA white dwarf cumulative age func¬ 
tion. The slope of a straight line fit to the last Gyr (red dashed 
line) gives the average DA WD formation rate. 


can be found in previous works jGarcf a-Berro et al.1 1 19991 : 
[Torres et al.ll2002 : iGarcfa-Berro et al.l 20041) . 

Any Monte Carlo code has at its very core the idea 
of repeated random sampling, i.e. generating the statistical 
properties of data from known distributions, that is used to 
obtain the initial characteristics (mass, time of birth, ini¬ 
tial position, kinematics, as well as other interesting quan¬ 
tities) of every star that will come to form part of the ini¬ 
tial synthetic population. One main ingredient here is the 
pseudo-rando m number ge nerator, for which we use the al¬ 
gorithm from ljamesi (Il990l ). This produces a uniform proba¬ 
bility density between (0,1) with a repetition period of over 
10^®, more than sufficient for most practical purposes. The 
next important step is using adequate probability distribu¬ 
tion functions for sampling the stellar properties. These dis¬ 
tribution functions are crucial inputs that define each par¬ 
ticular population. Thus, stellar m asses are samp led using 
the initial mass function (IMF) of iKrou 3 (l200ll) , a stan¬ 
dard choice, in partic ular considering the alleged universal 
character of the IMF dBastian et al.1 2010l) . The moment at 
which each star is born (tbom) is obtained in accordance 
with a star formation rate (SFR), assumed to be constant 
unless otherwise specified. The position of each star is ran¬ 
domly generated from a double exponential distribution of a 
constant Galactic scale height of 250 pc and a constant scale 
length of 3.5 kpc. The velocity distribution that we employ 
takes into account the differential rotation of the Galaxy, the 
peculiar veloci ty of the Sun and a sca le height dependent 
dispersion law ( Mihalas fc Binn^ll98ll ). Also, a metallicity 
value is assigned to each star acco rding to a Gaussian meta l- 
licity distribution as presented in iGasagrande et al.1 (1201 ll ). 

In order to reproduce the LSS-GAC, stars are only gen¬ 
erated in a cone delimited by —30° ^ 5 ^ 30° in Galactic 
latitude and 150° ^ I ^ 210° in Galactic longitude (Sec- 
tion[2]), with no restriction in terms of distance from the 
Sun. However, we define a test cone of up to 200 pc in length, 
in which we interactively examine the density of generated 
stellar mass until we re ach a limit density value f or the local 
stellar population (e.g. iHolmberg fc Flvnnl 120001 ). We scale 
this limit in order to obtain a final restricted WD sample of 
the same order as the observed one. 
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Table 2. The four models adopted in this work with the aim of 
reproducing the observed LSS-GAC DA WD population. Model 
1 is our standard model. 


Model 

SFR 

IFMR 

Slope for the 
massive regime 

1 

Constant 

Catalan et al. ^2008) 

0.10 

2 

Constant 

Catalan et al. f2008') 

0.06 

3 

Constant 

Ferrario et al. ('2005') 

0.10 

4 

Bimodal 

Catalan et al. f2008') 

0.10 


In a next step we set a 9.5 Gyr age for the thin disk 
(tdisk) and interpolate the main sequence lifetimes (Ims) of 
the generated stars using the BaSTI grids accord ing to stel¬ 
lar mass and metallicity (IPietrinferni et ani2004l . Knowing 
fdisk, fMS and fborn, we Can easily evaluate which of those 
stars have had time to become WDs. If that is the case the 
WD cooling age is simply given by tc = tdisk - Ims — tbom- 
Also, by knowing the mass of the WD progenitor we can 
compute the WD mass using an initial-to-final mass rela¬ 
tion (IFMR), which will be further detailed in Section lG^ 
WD luminosities, effective temperatures, surface gravities 
and UBVRI and Mboi magnitudes are then obtained by in¬ 
terpolating the inferred WD masses and cooling ages along 
the following cooling tracks: for WD masses smaller than 

1. 1 M[7) and larger tha n 0. 45 Mp we use the CO sequences 
of iRenedo et al] ll2010h and lAlthaus et ahl (l2010bll . while for 
W D masses abov e the u pper value we employ the O Ne tracks 
of lAlthaus et al.l ll2005l ') and lAlthaus et ^ (l2007l l. Finally, 
we convert the UBVR I magnitude s into the ugriz system 
using the equations of l.lordi et ahl (l200fih. taking into ac - 
count the 3D Galactic extinction map of [chen_et_alj (l2014l l 
and the extinction coefficients of lYuan et al.1 (l2013l l. 

Some of the DA WD stellar parameters derived from the 
fits to the LSS-GAC spectra are subject to relatively large 
errors (Section|4]). It is therefore necessary to account for 
those uncertainties in the simulated WD populations before 
comparing the synthetic and observational data sets. The 
effective temperature errors of the observed sample, which 
show a modest increase with increasing temperature, were 
fitted by a third order polynomial such that the error for rel¬ 
atively cool (~ 10, 000 K) WDs is about 300 K, increasing to 
~ 1, 000 K for WDs as hot as 25, 000 K. We adopt this poly¬ 
nomial relation for deriving effective temperature errors of 
our simulated WDs. The observational errors of log g clus¬ 
ter around ~ 0.2 dex, and we take this value as the surface 
gravity uncertainty of the synthetic WDs. The values of ef¬ 
fective temperature and surface gravity for each simulated 
WD are re-defined considering a random value within the er¬ 
ror range defined for the two quantities. We then interpolate 
new values of mass, luminosity, cooling age, and bolometric 
and absolute magnitudes from the redefined Teg and logg 
values. This results in, for example, an average error in mass 
of ~ 0.1 Mq for our simulated WDs. Photometric errors are 
also taken into account. They are directly derived from the 
photom etric uncertaint ies associated with the XSTPS-GAG 
survey llLiu et al.ll2014l L 


6.2 Models 

As presented in Section[5l the observational sample exhibits 
several specific features that are clearly visible in the LF, MF 
and cumulative AF of the LSS-GAG WD sample (Figs. [S] [9] 
an d[T0ll. We attempt to reproduce those features by employ¬ 
ing the above described population synthesis code. 

Given the apparent excess of massive WDs seen in the 
MF (Fig.[9)l, we attempt to reproduce this feature focusing 
on three parameters of the simulations that can affect the 
final WD mass distribution, namely the SFR, the IFMR 
up to an initial mass of about 6 Mq (the mass range of 
the zero-age main sequence for which WDs with CO cores 
are produced) and the slope of this (linear) relationship for 
the high-mass end (WDs with ONe cores). We start with 
our fiducial model, from now on called model 1, which uses 
a constant SFR, the piecewise linear IFMR introduced by 
ICatalan et al.1 ll2008ll for th e CO WD regime , and a slope of 
0.1 for the massive regime (llben et al.ll 19976 . 

We then consider three additional models, in which we 
vary only one of the above three parameters with respect to 
model 1. In mod el 2 we employ a sl ope of 0.06 for the IFMR 
of massive WDs (I Weidemannll200^') . The reason for lowering 
the slope of the relationship to this value is to expand the 
range of initial masses that can produce massive WDs (over 

1 Mq) in the hope of reproducing an excess. In order to en¬ 
sure the continuity of the IFMR over the entire WD mass 
range and to be consistent with the upper CO WD mass 
limit, we consider that all stars with masses between 6 and 
11 Mq become WDs of core masses ranging from 1.1 to 1.4 
Mq , which neatly gives this slope. Extending the mass range 
up to 1.4 Mq is probably wrong , given that a WD t hat mas¬ 
sive would most likely explode (iRitossa et al.lll999fl . but for 
the purposes of the current test is an a cceptable assumption . 
Model 3 uses the curved IFMR from iFerrario et al.1 (1200511 
for the CO regime, which, according to these authors, re¬ 
sults in a better agreement with the WD mass distribution 
as compared to when a linear fit is used. In model 4 we 
use the bimodal SFR of iRowell ll2013ll . which has two broad 
peaks at around 2 and 7 Gyr ago. This SFR should favor 
an increase in the number of massive WDs during the last 

2 Gyr given their shorter main sequence lifetimes. For each 
model we perform 10 individual realizations, and we com¬ 
pute the ensemble average of all the relevant quantities. A 
summary of the input parameters used for each model is 
given in Tabled 


6.3 The selection function 

Once the synthetic DA WD samples have been obtained 
for the different models outlined in the previous Section, it 
becomes necessary to evaluate which of those synthetic WDs 
would have been observed by the LSS-GAC survey. Here, we 
describe how the selection process is performed. 

The first step is to evaluate the effect of the LSS-GAC 
target selection criteria (Section[2]). This is done indepen¬ 
dently for each of the 10 realizations of each model. For 
this purpose, the g,r,i magnitudes of all WDs that are part 
of a given synthetic population are embedded within the 
XSTPS-GAG photometric catalogue and the selection crite¬ 
ria is then applied to the entire resulting population. The 
magnitude limits of the LSS-GAC survey are 14 mag ^ r ^ 
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Table 3. The synthetic WD populations are passed through con¬ 
secutive steps of filtering that gradually reduce the number of 
surviving objects. We show here an example for one of the 10 
performed realizations of our standard model (Table[2]l. In the 
last two columns we indicate the percentage of WDs that sur¬ 
vive respect to the previous step and the percentage of WDs that 
survive respect to the initial population, respectively. 



Filter 

Awd 

% 

% 


Initial sample 

3,874 




Initial sample within magnitude 
limits of the LSS-GAC 

2,132 



Step 1 

Selection Criteria 

1,748 

82 

82 

Step 2 

GAC plates -h fibre allocation 

192 

11 

9 

Step 3 

S/N ^5 

102 

53 

5 

step 4 

Completeness -t- spectral fit 

77 

75 

3.5 


18.5mag (Section[^, and within those limits the LSS-GAC 
criteria efficiently selects 80-85 per cent of the total simu¬ 
lated WD population, depending on the model. This frac¬ 
tion increases to 85-90 per cent if we consider 14 mag ^ r ^ 
18 mag. The high success rate of selecting WDs is not unex¬ 
pected, considering that the LSS-GAC survey is specifically 
developed to efficiently target stars of all colours, including 
blue objects such as WDs (Section[2I. 

In a second step we evaluate which of the simulated 
WDs that are selected by the LSS-GAC criteria fall within 
the field of view (5deg^) of the LSS-GAC plates actually ex¬ 
ecuted (Table Al). If this is the case, an additional condition 
is that the simulated WDs are required to fulfill the mag¬ 
nitude limits of the plates/spectrographs, otherwise they 
would not have been observed. In practice, we consider the 
distances between the position defined by the right ascen¬ 
sion and declination of each simulated WD and the central 
positions of the 16 spectrographs of the plate where the syn¬ 
thetic WD falls (also defined by their right ascensions and 
declinations) and evaluate whether or not the r magnitude 
of the synthetic WD is within the magnitude limits of the 
nearest spectrograph. If all those conditions are fulfilled, we 
then consider the probability of a given target to be allo¬ 
cated a fibre (some fibres are used for sky observations). 
This probability is simply given by Wpec/(A^spec + A^sky), 
and is generally ~0.9. Wpec is the number of target spectra 
observed by the spectrograph, and Wky is the number of 
fibres allocated for sky observations. 

If the synthetic WDs survive all the previously ex¬ 
plained filtering process we consider the LSS-GAC survey 
would have observed them. Therefore, in a third step we 
consider the probability for each simulated WD to have a 
LAMOST spectrum of S/N ratio > 5 in both the blue and 
red arms. For each synthetic WD we calculate the fluxes 
from their associated g and r magnitudes, add and subtract 
a 5 per cent of flux in each case and calculate the mag¬ 
nitudes that result from this exercise {g+, g- and r+, r_; 
where the sufixes -I- and — indicate that we have added and 
subtracted the 5 per cent of the corresponding flux). We then 
consider all targets observed by the respective spectrograph 
(i.e. the spectrograph where the simulated WD falls) having 
g- < g < g+ and r_ < r < r+, and calculate the median 
S/N ratio of their LSS-GAC spectra in the two bands. If no 


observed spectra are found satisfying the above magnitude 
ranges, or if one of the median S/N ratios is smaller than 5, 
the synthetic WD is then excluded from the analysis. This 
exercise takes into account nigh-to-night variations of S/N 
ratio that may arise e.g. from varying observing conditions, 
as the S/N is evaluated specifically for objects observed dur¬ 
ing the same night with the same plate/spectrograph, and 
of similar magnitudes as the simulated WD of concern. 

Finally, in a fourth step we take into account the spec¬ 
troscopic completeness of the observed sample (the fraction 
of LSS-GAC DA WDs that we have identified among all 
DA WDs observed) as well as consider the fact that we are 
not able to obtain reliable stellar parameters for 19 per cent 
of the observed sample. We have estimated a spectroscopic 
completeness of 95 per cent ('Section l3.3ll . Therefore, we ran¬ 
domly exclude 5 per cent of all synthetic WDs that passed 
the previous filters. After this correction, we proceed by ran¬ 
domly excluding 19 per cent of the surviving systems. 

In order to minimize the effects of the random exclusion 
of synthetic WDs, we repeat steps two to four 20 times per 
model realization. Given that each of the four models consid¬ 
ered (Table[2)) counts 10 realizations, we obtain 200 different 
final synthetic populations for each model. The number of 
simulated WDs that pass the entire selection process de¬ 
scribed above vary slightly from model to model (and real¬ 
ization to realization) and yields synthetic samples of 65-85 
objects, similar to the number of WDs in the observed sam¬ 
ple, 75 DA WDs. An example of how the number of synthetic 
WDs gradually decreases as they are passed through each of 
the filters of our selection process is shown in Table[31 In a 
final step we use bootstrapping techniques to produce syn¬ 
thetic samples of the same number of objects as the observed 
one. 

The final LF, MF and cumulative AF for each model 
are the result of averaging 200 individual functions derived 
from each of the independent realizations. These are shown 
in Figs. llllll2l and n^ respectively, where we also include the 
LF, MF and cumulative AF derived from the observational 
sample. A comparison of the synthetic and observed func¬ 
tions is presented and discussed in detail in the following 
Section. 


7 DISCUSSION 

In this section we compare the LFs, MFs and cumulative 
AFs (i.e. birth rates) that result from our numerical simula¬ 
tions to those derived observationally. Before comparing in 
detail the simulated and the observed distributions, we com¬ 
pare the WD populations obtained from each of the models 
employed here. We also estimate the number of DA WDs 
that the LSS-GAC will eventually observe. 

7.1 The final expected number of LSS-GAC DA 
white dwarfs 

The LSS-GAC selection criteria (Section[21) applied to our 
simulated WD populations results in ~80-85 per cent of the 
synthetic DA WDs falling within the magnitude limits of 
the survey being selected for observations fSection l6.3D . The 
total number of WDs generated by each model oscilates be¬ 
tween 3,800 and 3,900, which reduces to 2,130-2,160 if we 
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Figure 11. Simulated DA WD LFs (grey solid lines) resulting from the models considered in Section|6]2] (see Table l^ . The observed 
LSS-GAC WD LF is shown as a black solid line. 


consider the magnitude limits of the LSS-GAC. This implies 
that, on average, >^1,700-1,850 DA WDs could be potentially 
observed at the end of the survey, depending on the assumed 
model. 


We have shown that ~50 per cent of our synthetic 
DA WDs would have S/N^5 if observed by the LSS-GAC 
(see TablejS]). This percentage is expected to increase up 
to ~2/3 for the data release 2 (and further releases) of 
LSS-GAC spectra (private communication). The number 
of LSS-GAC DA white dwarfs expected to have spectra 
of S/N^5 at the end of the survey is thus Ada =(1,700- 
1,850 — Nohs) X 2/3 ~ 1,070-1,170, where Aobs = 92, i.e. 
the number of currently observed DA WDs. Considering 
spectroscopic completeness and spectral fitting effects (Sec- 
tion l6.3l Table[3ll, which exclude ~ 25 per cent of the DA 
WD spectra with S/N^5, the final number of LSS-GAC DA 
white dwarfs with available and reliable stellar parameters 
at the end of the survey is expected to be ~800-875, i.e. 
approximately one order of magnitude higher than the cur¬ 
rent number of observed LSS-GAC DA white dwarfs with 
reliable stellar parameters. 


7.2 Effects of observational uncertainties 

We have employed four different models to simulate the WD 
population in the Galactic anti-center with the aim of con¬ 
straining what set of assumptions (SFR, IFMR,...) fits bet¬ 
ter the observational data. As expected, the intrinsic prop¬ 
erties of the simulated populations differed from model to 
model. However, these properties are altered when the ob¬ 
servational uncertainties are incorporated (Section l 6 .ll) . For 
example, the simulated mass distributions become broader 
and lose detail, and more importantly, peak at larger values 
(e.g. the median of the distribution shifts from 0.57 Mq to 
0.60 M 0 for Model 2). Hence, the incorporation of obser¬ 
vational uncertainties results in less prominent differences 
between the synthetic WD parameter distributions. 

This effect is enhanced when we take into account the 
selection biases. In order to illustrate these effects together, 
we show in Fig.Qthe correlations between the effective tem¬ 
peratures, masses, cooling ages and distances for the syn¬ 
thetic population (red stars) and compare them to those 
obtained from the observational sample (black dots). For 
the seek of clarity we chose one typical realization of our 
Model 1, although very similar results are obtained for the 
other realizations and models. It becomes obvious that the 
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Figure 12. Same as Fig. llll but for the MFs. The fraction of massive WDs (> 0.8 Mq) for both the observed and the simulated samples 
are provided in the top right of each panel. 


model reproduces well the observational data, and that the 
correlations between the considered parameters follow the 
same pattern as the observational one (Section|4]). 


7.3 The luminosity function 

Fig.[TT]shows the LFs derived from our simulated samples, as 
well as that deduced from the observational data. The uncer¬ 
tainties in the simulated functions were derived in the same 
way as for the observed one fSection l5.ll) . The space density 
obtained for models 1, 2, 3 and 4, are 0.96±0.19, 0.98±0.21, 
1.16 ± 0.20 and 1.06 ± 0.20 x 10“® pc“® respectively. Al¬ 
though these values are slightly higher than the space den¬ 
sity derived from our observations (0.83 ±0.16 x 10“® pc~^, 
Section[3, they perfectly match within the error bars. It is 
evident that there is an overall good agreement (within the 
error bars) between the simulated and the observed LFs, 
except at Mboi 2 and 6 mag, where the observed LF pre¬ 
dicts a considerably higher space density. It has to be noted, 
however, that the number of targets falling within bins of 
Afboi < 7 mag is small (18 per cent of the total observed 
sample). Hence, the observed LF in those high luminosity 
bins is subject to low number statistics and the apparent 


increase of the observed LF at those specific bins should 
be taken with some caution. A further inspection of Fig. 1111 
reveals that, because of the reasons explained above (Sec- 
tion l7.2l) . no model seems to have an obvious advantage in 
reproducing the observational data. 


7.4 The mass function 

After applying the LSS-GAC target selection criteria and 
the target selection process described in Section [6l3l to our 
simulated populations, the MFs yielded by all simulations 
are rather similar fFig. ll2ll . In addition, synthetic (single) 
WDs of masses as low as 0.35 Mq are now possible as a con¬ 
sequence of incorporating observational uncertainties. This 
effect partly explains the apparent over-abundance of low- 
mass (< 0.45 Mq) WDs in the observed MF (black solid line 
in Fig. mi). Alternatively, a relative large fraction of low- 
mass WDs in the ob served sample could be the res ult of 
binary star evolution dRebassa-Mansergas et aT]|201ll l. The 
companions are likely to be cooler and/or more massive 
WDs, or low-mass late-type main sequence stars, although 
other exotic companions such as brown dwarfs cannot be 
ruled out. 
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Figure 13. Normalised LSS-GAC DA white dwarf MF (gray) 
and normalised mass dis tribution of DA WDs from the local and 
volume-limited sample of iGiammichele et alj ll2012t) (red dashed). 


It is also clear that none of onr models manages to com¬ 
pletely reproduce the observed behaviour at high mass bins, 
i.e. the fraction of massive WDs (> 0.8 Mq) relative to those 
of typical mass (0.6 Mq) is higher in the observed sample. 
We discuss possible scenarios leading to this feature below. 


7.4-1 The initial-to-final mass relation 

The currently available IFMRs have been derived from ob¬ 
servational data that exhibit large scatter in the initial-to- 
final mass diagram (see for example Fig. 1 of I Catalan et al.l 
l2008h . In one of our models we have investigated the ef¬ 
fect of varying the slope of the IFMR for producing a wider 
range of massive WDs (see Section[S2] and Table[5]). We 
have also explored t he effect of employing a curved IFMR 
llFerrario et al.ll2005l) . To further investigate the impact of 
the large scatter in the initial-to-final mass diagram to the 
simulated MF, two additional models (models 5 and 6) are 
d eveloped that t a ke int o account the error bars of the IFMR 
of I Catalan et al.l ll2008l ) so that the IFMR is virtually moved 
“up” in one model and “down” in the other model. The re¬ 
maining free parameters of models 5 and 6 are the same as 
for our standard model (Tabled . The results show that the 
MFs obtained from these two models do not differ signifi¬ 
cantly from those shown in Fig.[T2]and therefore are not able 
of reprodncing the high-mass excess present in our observed 
MF. 


7 . 4.2 S/N ratio and 3D model atmosphere correction 
effects 

Two additional plansible explanations for the large fraction 
of massive WDs observed are effects of limited S/N ratios 
and 3D model atmosphere corrections. 

The LSS-GAC spectra considered in this work have a 
minimum S/N ratio of 5. It is therefore possible that some 
systematic uncertainties in the WD stellar parameters result 
as a consequence of the relatively low S/N ratio of some WD 
spectra. This may lead to the masses of some WDs being 


overestimated. In order to investigate this possibility we re¬ 
derive the observed MF excluding all systems with spectra 
of a S/N ratio below 8. This leaves us with 52 DA WDs. We 
decided not to increase the S/N threshold to higher values 
because otherwise the number of massive (fainter and with 
systematically lower S/N ratios) WDs that would survive 
the cut would be severely reduced. The MF that results from 
this exercise does not differ significant from the one obtained 
using the full sample, and displays as well a large fraction 
of massive WDs. Therefore, S/N ratio effects are unlikely to 
be the cause of the excess of massive WDs observed. 

The DA WD sample analysed in this work includes 
cool WDs for which we have applied the 3D model atmo¬ 
sphere corrections to their stellar parameters deduced from 
ID model atmosphere fitting: effective temperature, surface 
gravity, and hence mass. If those corrections are somehow 
incorrect, they may lead to an apparent overabundance of 
massive WDs. To explore this possibility we re-derived the 
MF excluding all WDs in our sample with an effective tem¬ 
perature below 13,000 K. This results in a sub-sample of 
58 DA WDs. The MF deduced from this sub-sample again 
presents a clear overabundance of massive WDs. We there¬ 
fore find that the overabundance of massive WDs is unlikely 
caused by the possible effects related to the 3D model at¬ 
mosphere corrections. 


7 . 4.3 Effective temperature and surface gravity error 
effects 

The effective temperature and surface gravity errors ob¬ 
tained fitting the Balmer lines are not independent. The 
reason for this correlation is that the strength of the Balmer 
lines is largely determined by the ionization balance. That 
is, if the equivalent width of e.g. H/3 is fixed, a higher as¬ 
sumed effective temperature will need a higher surface grav¬ 
ity, as the higher pressure is needed to compensate. The line 
shape is the second order effect, which determines where on 
the correlation line the best solution is. In order to ives- 
tigate whether or not this effect may explain the excess of 
high-mass WDs observed, our simulations should have taken 
into account not only the errors in these quantites (see Sec- 
tions l6.3l and l7.2ll . but also their correlation. 

However, the strength of this correlation not only de¬ 
pends on the effective temperature range, but also the first 
order change of line strength disappears, and with it the 
correlation, near the maximum strength of the line. More¬ 
over, the correlation is not always apparent for small errors. 
Hence, quantifying the correlation between the effective tem¬ 
perature and surface gravity errors is a notable endeaveour, 
which is beyond the scope of this paper. Thus, whether or 
not such a correlation may explain the apparent excess of 
massive WDs remains an open question. 

7 . 4.4 WD+WD mergers 

An exciting possible explanation for the excess of high- 
mass WDs in the observed MF is that a relatively large 
fraction of tho se stars are the result of merg ers of two 
low-mass WDs llMarsh et al.lll99^ : IVennesiri999i) . Although 
no population synthesis study hitherto predicts such a 
large fraction of high-mass WDs as the outcome of WD 
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Figure 14. Same as Fig. llll but for the cumulative AFs. 


mergers (e.g. Han et al.iri994l : lHarJll99^ : iMeng et al.ll2008l : 
iToonen et al.l 20121 '). this scenario has been ad opted in some 
of the previous observational studies (e.g. iLiebert et al.l 
1200111 : iGiammichele et al.ll2012l 'l. To further investigate this 
hypothesis we compare in Fig.[T3] the normalised MF ob¬ 
tained in this work with the normalised mass distribu- 
tion of DA WDs i n the local, volume-limited sample of 
IGiammichele et al.l (l2012f ). In becomes obvious that the peak 
at 0.6 Mq is less pronounced in our normalised MF, an ef¬ 
fect likely related to the fact that we are subject to larger 
observational uncertainties which broaden the distribution. 
Interestingly, whilst the high-mass peak in th e normalised 
mass distribution of IGiammichele et al.l l|2012l ~l is found at 
1 Mq, our normalised MF shows two apparent peaks at the 
0.9 and 1.1 M© bins and reflects a scarcity of systems at 
the 1 Mq bin. Although this discrepancy is likely due to 
our larger uncertainties, which are sufficient to shift objects 
across bins, both studies favour the hypothesis that an ex¬ 
cess of massive 'WDs seems to exist. 

If this excess of massive WDs may arise as a conse¬ 
quence of WD binary mergers, then the high merger rate 
implied by the observed excess of massive WDs may indicate 
a much more important role of the double-de generate chan- 
nel for the production of Type la supernovae JWang fc HanI 


I2OI2I : iToonen et al.l[2ni2h . Analysing the expected merger 
rates of WDs based on our observations and comparing them 
to the observed rates of Type la supernovae thus seems to 
be a worthwhile exercise, and we will pursue this elsewhere. 
High-mass W Ds that result from merg ers are expected to 
be magnetic llGarcfa-Berro et al.ll201^ . therefore we may 
expect to find signs of magnetic fields in our massive WDs 
that should help test this hypothesis. 

The large WD merger rate suggested by the current 
work indicates that an even larger number of close WD bi¬ 
naries may exist in the Galaxy that have not yet merged. 
Those close binaries could be a main source of gravitational 
waves to be detected by future faciliti es such as the space 
interferometer eLISA dNelemansI 1201,11 1 . Therefore, indirect 
support in favour of the merger scenario may be obtained by 
analysing the population of close double WDs that eLISA 
will discover. 


7.5 The average DA WD formation rate 

The cumulative AFs derived from our simulated populations 
are illustrated in Fig. lldi where the observed cumulative AF 
is also displayed. There is an overall good agreement be¬ 
tween our simulations and the osbervations for coolinjg ages 
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up to 1 Gyr, except perhaps for our model 4, which seems 
to systematically overestimate the space density for cooling 
age bins <0.4 Gyr (note that in model 4 we are employing 
a bimodal star formation rate). For cooling ages larger than 
1 Gyr the discrepancies between the models and the obser¬ 
vations arise due to the scarcity of WDs at those specific 
cooling ages. 

Fitting the simnlated cumulative AFs with a straight 
line (see Section [531l we derive average DA WD forma¬ 
tion rates of 6.04±0.05, 6.42±0.05, 5.85±0.02 and 5.97±0.04 
xl0“^® pc~® for models 1, 2, 3 and 4 respectively. 

These values agree with the average formation rate derived 
from our observations within the errors (5.42 ±0.08 x 10“^® 
pc“® yr~^, Section lS^ . Because of the observational uncer¬ 
tainties fSection l7.2ll . we find that no model seems to have 
an obvious advantage in reproducing the observational data. 


8 SUMMARY AND CONCLUSIONS 

The recently initiated LAMOST Spectroscopic Survey of the 
Galactic anti-center, the LSS-GAC, selects targets for spec¬ 
troscopic observations following a well-defined criteria. This 
significant advantage over previous surveys has allowed us to 
present a well-characterised magnitude-limited sample of 92 
LSS-GAG hydrogen-rich (DA) white dwarfs from the data 
release 1. Our catalogue is expected to be ~ 95 per cent com¬ 
plete. We have determined the stellar parameters (surface 
gravity, effective temperature and mass), absolute and bolo- 
metric magnitudes, and distances of 75 DA white dwarfs. 
Taking into account volume corrections we have derived an 
absolute lower limit for the space density of DA white dwarfs 
of 0.83T0.16 X 10“® pc“^. We have also obtained prelimi¬ 
nary observed LSS-GAC DA white dwarf luminosity, mass 
and cumulative age functions. The luminosity function re¬ 
sembles those found in previous observational studies. The 
mass function reveals an excess of massive white dwarfs. Fi¬ 
nally, the DA white dwarf formation rate derived from the 
cumulative age function is 5.42±0.08x pc~® in 

good agreement with other recent studies. 

We have simulated the DA white dwarf population in 
the Galactic anti-center using an existing Monte Carlo code 
adapted to the characteristics of the LSS-GAC. For this pur¬ 
pose, and specially aiming at reproducing the observed ex¬ 
cess of massive white dwarfs, we have employed four differ¬ 
ent models. All those models take into account the obser¬ 
vational uncertainties, both spectroscopic (i.e., we incorpo¬ 
rate errors in the stellar parameters of our simulated white 
dwarfs based on the observational errors) and photometric. 
We find that the LSS-GAC criteria selects >^80-85 per cent 
of all simulated white dwarfs with 14 sC r ^ 18.5 mag (the 
magnitude limits of the survey) in each model, thus provid¬ 
ing robust evidence for the high efficiency of LSS-GAC in 
targeting white dwarfs. Once the observational uncertainties 
have been taken into account in our simulations, the distri¬ 
bution of stellar parameters are similar for all models. We 
find that all our simulations reproduce well the observed lu¬ 
minosity function, however no particular model seems to fit 
better the data. 

None of our considered models is able to reproduce the 
observed excess of massive DA white dwarfs. We have inves¬ 
tigated possible explanations for this feature and concluded 


that a plausible scenario is that a sizable fraction of those 
massive white dwarfs are products of mergers of two initially 
lower-mass white dwarfs. If that is the case, then the white 
dwarf merger rate in our Galaxy is considerably higher than 
currently assumed. This may have important implications 
for the production of Type la supernovae via the double¬ 
degenerate channel. 

Finally, it is important as well to emphasise that al¬ 
though our study represents an important step forward 
towards unveiling the underlying population of DA white 
dwarfs in the Galaxy, the size of the LSS-GAC sample is 
small, and that the stellar parameters we derived for some 
objects are subject to relatively large uncertainties. Forth¬ 
coming LSS-GAC data releases are expected to increase the 
number of DA white dwarfs by one order of magnitude. In 
addition, the quality of the LAMOST spectra will improve, 
which will reduce the uncertainties in the stellar parame¬ 
ter determinations. We will hence derive updated luminosity 
and mass functions and DA white dwarf formation rates at 
the end of the survey. 
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Table Al. LSS-GAC plates, including associated right ascensions and declinations, used for observations. The date of observations are 
also included. 


Date 

Plate 

RA 

[deg] 

DEC 

[deg] 

Date 

Plate 

RA 

[deg] 

DEC 

[deg] 

2011-10-03 

PA09B_kedal 

63.79331 

29.90206 

2012-10-12 

GAC090N33B1 

90.06008 

33.13693 

2011-10-03 

PA09B_keda2 

63.79331 

29.90206 

2012-10-13 

GAC072N32B1 

72.32949 

32.58819 

2011-10-03 

PA09M_kedal 

63.79331 

29.90206 

2012-10-13 

GAC072N32M1 

72.32949 

32.58819 

2011-10-03 

PA09M_keda2 

63.79331 

29.90206 

2012-10-17 

GAC083N27B1 

83.98125 

27.66235 

2011-10-05 

PA09F_kedal 

63.79331 

29.90206 

2012-10-17 

GAC083N27B2 

83.98125 

27.66235 

2011-10-05 

PA09M_kedal 

63.79331 

29.90206 

2012-10-18 

GAC085N33B2 

85.74273 

33.31455 

2011-10-05 

PA09M_keda2 

63.79331 

29.90206 

2012-10-19 

GAC080N32B1 

80.00299 

32.78560 

2011-10-21 

PAIB 

45.87566 

28.26991 

2012-10-19 

GAC080N32F1 

80.00299 

32.78560 

2011-10-21 

PB8B 

91.59353 

29.51241 

2012-10-19 

GAC080N32M1 

80.00299 

32.78560 

2011-10-21 

PB8M 

91.59353 

29.51241 

2012-10-24 

GAC051N24B1 

51.07698 

24.72406 

2011-10-28 

GAC.060N28.B1 

60.35798 

28.50189 

2012-10-24 

GAC051N24M1 

51.07698 

24.72406 

2011-10-28 

GAC.105N29.B1 

105.02926 

29.77213 

2012-10-24 

GAC100N32B1 

100.68762 

32.55887 

2011-11-08 

GAC.097N28.B1 

97.59171 

28.21211 

2012-10-24 

GAC100N32M1 

100.68762 

32.55887 

2011-11-08 

GAC.122N29.B1 

122.39659 

29.09308 

2012-10-25 

GAC086N24B1 

86.98389 

24.68684 

2011-11-09 

GAC.072N28.B1 

72.31558 

28.35083 

2012-10-25 

GAC086N24B2 

86.98389 

24.68684 

2011-11-09 

GAG.101N28.B1 

101.18941 

28.97093 

2012-10-25 

GAC086N24M1 

86.98389 

24.68684 

2011-11-09 

GAC.118N28.B1 

118.33777 

28.05199 

2012-10-27 

GAC090N26B1 

90.89135 

26.52913 

2011-11-10 

GAC.089N28.B1 

89.14070 

28.94227 

2012-10-27 

GAC090N26B2 

90.89135 

26.52913 

2011-11-10 

GAG.089N28.B2 

89.14070 

28.94227 

2012-10-27 

GAC102N27B1 

102.29246 

27.19037 

2011-11-10 

GAC.089N28.B3 

89.14070 

28.94227 

2012-10-27 

GAC113N27B1 

113.05368 

27.12516 

2011-11-10 

GAC.113N28.B1 

113.63163 

28.68658 

2012-10-29 

GAC078N26B1 

78.08962 

26.45461 

2011-11-11 

GAC.082N29.B1 

82.41940 

29.18646 

2012-10-29 

GAC098N33B1 

98.42774 

33.02404 

2011-11-12 

GAG.080N28.B1 

80.84530 

28.93676 

2012-10-29 

GAC113N26B1 

113.98061 

26.89574 

2011-11-14 

GAC.087N27.B1 

87.66559 

27.50503 

2012-10-29 

GAC117N24B1 

117.32855 

24.48873 

2011-11-14 

GAC.106N28.B1 

106.85366 

28.17669 

2012-10-31 

GAC097N26B1 

97.23493 

26.96746 

2011-11-20 

GAC.063022N281243.fi 

97.59171 

28.21211 

2012-10-31 

GAC114N33B1 

114.58637 

33.18676 

2011-11-20 

GAC.063022N281243.M1 

97.59171 

28.21211 

2012-11-06 

test_055N28.Bl 

55.32663 

28.70276 

2011-11-20 

GAC.080935N290534.M1 

122.39659 

29.09308 

2012-11-06 

test_080N33.Bl 

79.75012 

33.74839 

2011-11-23 

GAC.107N27.B1 

108.98819 

27.89742 

2012-11-06 

test_080N33.B2 

79.75012 

33.74839 

2011-11-23 

GAC.107N27.M1 

108.98819 

27.89742 

2012-11-06 

test_122N25.Bl 

122.55461 

25.84448 

2011-11-24 

GAC.067N28.M1 

72.31558 

28.35083 

2012-11-07 

GAC049N32B1 

49.14662 

32.18402 

2011-11-24 

GAC.106N28.M1 

106.85366 

28.17669 

2012-11-07 

GAC049N32M1 

49.14662 

32.18402 

2011-11-26 

GAC.082N27.M1 

87.66559 

27.50503 

2012-11-07 

GAC069N25B1 

69.84645 

25.21827 

2011-11-29 

GAC.096N28.M1 

100.33707 

28.19664 

2012-11-07 

GAC107N32B1 

107.68018 

32.61873 

2011-11-30 

PBOIB 

47.40309 

29.07708 

2012-11-07 

GAC117N27B1 

117.47942 

27.36318 

2011-12-03 

GAC.045N28.B1 

45.87566 

28.26991 

2012-11-13 

GAC073N32B1 

73.87681 

32.78394 

2011-12-03 

GAC.045N28.M1 

45.87566 

28.26991 

2012-11-13 

GAC073N32M1 

73.87681 

32.78394 

2011-12-03 

GAC.083N28.M1 

89.14070 

28.94227 

2012-11-13 

GAC102N32B1 

102.42213 

32.60676 

2011-12-07 

GAC.106N28.B1 

106.85366 

28.17669 

2012-11-13 

GAC102N32M1 

102.42213 

32.60676 

2011-12-11 

GAC.045N28.B1 

45.87566 

28.26991 

2012-11-14 

GAC098N33F1 

98.42774 

33.02404 

2011-12-11 

GAC.065N28.B1 

68.65830 

28.96115 

2012-11-14 

GAC098N33M1 

98.42774 

33.02404 

2011-12-11 

GAC.079N29.B1 

82.52563 

29.54832 

2012-11-14 

GAC121N33B1 

121.03494 

33.03074 

2011-12-14 

GAC.078N28.B1 

80.84530 

28.93676 

2012-11-17 

GAC079N24B1 

79.38239 

24.00968 

2011-12-14 

GAC.078N28.M1 

80.84530 

28.93676 

2012-11-17 

GAC079N24M1 

79.38239 

24.00968 

2011-12-14 

GAC.105N29.B1 

105.02926 

29.77213 

2012-11-19 

GAC058N31F1 

58.01852 

31.16858 

2011-12-15 

GAC.082N27.B1 

87.66559 

27.50503 

2012-11-19 

GAC121N33F1 

121.03494 

33.03074 

2011-12-16 

PB03B 

56.15493 

27.89745 

2012-11-19 

GAC121N33M1 

121.03494 

33.03074 

2011-12-16 

PB03M 

56.15493 

27.89745 

2012-11-22 

GAC049N27B1 

49.98249 

27.07113 

2011-12-17 

GAC.034118N284209.fi 

55.32663 

28.70276 

2012-11-22 

GAC089N24F1 

89.23383 

24.24969 

2011-12-18 

GAC.04h29.Bl 

61.75189 

29.00130 

2012-11-22 

GAC089N24M1 

89.23383 

24.24969 

2011-12-18 

GAC.04h29.Ml 

61.75189 

29.00130 

2012-11-23 

GAC083N24B1 

83.37816 

24.62881 

2011-12-19 

GAC.079N29.M1 

82.41939 

29.18646 

2012-11-23 

GAC104N26B1 

104.69754 

26.08108 

2011-12-20 

GAC.067N28.fi 

72.31558 

28.35083 

2012-11-23 

GAC120N25B1 

120.23281 

25.39284 

2011-12-21 

GAC.107N27.B1 

108.98819 

27.89742 

2012-11-25 

GAC117N31B1 

117.75082 

31.61353 

2011-12-22 

GAC.082N29.M1 

82.52563 

29.54832 

2012-12-04 

GAC041N29B1 

41.91085 

29.67853 

2011-12-23 

GAC.118N28.fi 

118.33777 

28.05199 

2012-12-04 

GAC117N24M1 

117.32855 

24.48873 

2011-12-23 

GAC.118N28.M1 

118.33777 

28.05199 

2012-12-05 

GAC065N31M1 

65.04145 

31.95317 

2011-12-24 

GAC.068N28.fi 

68.65830 

28.96115 

2012-12-06 

test_091N23.Bl 

91.70276 

23.63860 

2011-12-24 

GAC.068N28.M1 

68.65830 

28.96115 

2012-12-06 

test_091N23.Ml 

91.70276 

23.63860 

2011-12-24 

GAC.106N28.B1 

106.85366 

28.17669 

2012-12-06 

test_lllN36.Bl 

111.07162 

36.31091 

2011-12-25 

GAC.089N28.B3 

89.14070 

28.94227 

2012-12-06 

test_128N24_Bl 

127.87716 

24.08111 

2011-12-25 

GAC.089N28.M3 

89.14070 

28.94227 

2012-12-06 

test_128N24_B2 

127.87716 

24.08111 

2011-12-25 

GAC.113N28.M1 

113.63163 

28.68658 

2012-12-08 

GAC074N27B1 

74.97396 

27.32560 
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Table Al. Continued. 


Date 

Plate 

RA 

[deg] 

DEC 

[deg] 

Date 

Plate 

RA 

[deg] 

DEC 

[deg] 

2011-12-26 

GAC.047N29.B1 

47.40309 

29.07708 

2012-12-08 

GAC074N27M1 

74.97396 

27.32560 

2011-12-26 

GAC.091N29.B1 

91.59353 

29.51241 

2012-12-09 

GAC084N26B1 

84.73905 

26.61807 

2011-12-26 

GAG.091N29.M1 

91.59353 

29.51241 

2012-12-09 

GAC084N26M1 

84.73905 

26.61807 

2011-12-27 

GAG.105N29.M1 

105.02926 

29.77213 

2012-12-09 

GAC114N32B1 

114.97536 

32.00973 

2011-12-28 

GAG.108N27.M1 

108.98819 

27.89742 

2012-12-21 

GAC120N25M1 

120.23281 

25.39284 

2011-12-31 

PA06F 

68.65830 

28.96115 

2012-12-22 

GAC108N24B1 

108.10991 

24.12859 

2011-12-31 

PA06M 

68.65830 

28.96115 

2013-01-04 

test.076N27.Bl 

76.15819 

27.69602 

2012-01-02 

PB06B 

82.41939 

29.18646 

2013-01-04 

test_076N27.Ml 

76.15819 

27.69602 

2012-01-02 

PB06M 

82.41939 

29.18646 

2013-01-04 

test.094N28.Bl 

94.99571 

28.42676 

2012-01-03 

GAG.068N28.B1 

68.65830 

28.96115 

2013-01-04 

test_126N31.Bl 

126.02454 

31.30104 

2012-01-03 

GAC.101N28_M1 

100.33707 

28.19664 

2013-01-05 

GAC087N32B1 

87.85729 

32.12469 

2012-01-04 

GAG.063N29.B1 

63.79331 

29.90206 

2013-01-05 

GAC087N32M1 

87.85729 

32.12469 

2012-01-04 

GAC.106N28.M1 

106.85366 

28.17669 

2013-01-06 

GAC096N32B1 

96.30344 

32.27176 

2012-01-11 

GAG.050N29.B1 

50.08484 

29.04846 

2013-01-06 

GAC096N32M1 

96.30344 

32.27176 

2012-01-12 

GAG.045N28.B1 

45.87566 

28.26991 

2013-01-07 

GAC085N33M1 

85.74273 

33.31455 

2012-01-12 

GAG.100N28.B1 

100.33707 

28.19664 

2013-01-07 

GAC122N25B1 

122.55461 

25.84448 

2012-01-13 

GAG.061N29.B1 

61.75189 

29.00130 

2013-01-07 

GAC122N25M1 

122.55461 

25.84448 

2012-01-14 

PA09B1 

89.14070 

28.94227 

2013-01-08 

GAC086N24M2 

86.98389 

24.68684 

2012-01-14 

PA09M1 

89.14070 

28.94227 

2013-01-09 

GAC085N31B1 

85.14961 

31.35820 

2012-01-14 

PB02B 

52.19438 

30.37534 

2013-01-09 

GAC085N31M1 

85.14961 

31.35820 

2012-01-15 

PAOIB 

45.87566 

28.26991 

2013-01-11 

GAC094N27M1 

94.58662 

27.21015 

2012-01-15 

PA09B1 

89.14070 

28.94227 

2013-01-12 

GAC054N25M1 

54.29599 

25.99110 

2012-01-15 

PB05B 

70.09451 

29.97223 

2013-01-12 

GAC078N26M1 

78.08962 

26.45461 

2012-01-21 

GAC.060N28.fi 

60.35798 

28.50189 

2013-01-12 

GAC109N30M1 

109.51705 

30.95587 

2012-01-21 

GAC.100N28.M1 

100.33707 

28.19664 

2013-01-13 

GAC080N33B2 

79.75012 

33.74839 

2012-01-22 

GAC.063N29.M1 

63.79331 

29.90206 

2013-01-13 

GAC080N33M1 

79.75012 

33.74839 

2012-01-22 

GAC.089N28.F1 

89.14070 

28.94227 

2013-01-16 

GAC045N26B1 

45.47559 

26.46235 

2012-01-23 

GAC.061N29.M1 

61.75189 

29.00130 

2013-01-16 

GAC055N32B1 

55.61516 

32.93923 

2012-01-23 

GAC.089N28.F2 

89.14070 

28.94227 

2013-01-16 

GAC065N31B1 

65.04145 

31.95317 

2012-01-24 

GAC.060N28.M1 

60.35798 

28.50189 

2013-01-17 

GAC040N27B1 

40.86297 

27.70715 

2012-01-24 

GAC.080N28.M1 

80.84530 

28.93676 

2013-01-17 

GAC105N24B1 

105.60325 

24.21545 

2012-01-24 

GAC.108N27.M1 

108.98819 

27.89742 

2013-01-29 

GAC058N25B1 

57.85539 

25.16293 

2012-01-25 

GAC.052N30.M1 

52.19438 

30.37534 

2013-02-01 

GAC076N33B1 

76.53483 

33.91869 

2012-01-25 

GAC.097N28.fi 

97.59171 

28.21211 

2013-02-04 

GAC046N25B1 

46.36121 

25.25517 

2012-01-26 

GAC.055N28.M1 

55.32663 

28.70276 

2013-02-07 

GAC054N25B1 

54.29599 

25.99110 

2012-01-26 

GAC.089N28.M2 

89.14070 

28.94227 

2013-02-08 

GAC046N25B1 

46.36121 

25.25517 

2012-01-26 

GAC.122N29.M1 

122.39659 

29.09308 

2013-02-08 

GAC062N26B1 

62.70775 

26.48095 

2012-01-29 

PAOIM 

45.87566 

28.26991 

2013-02-08 

GAC109N31B1 

109.37725 

31.69805 

2012-01-29 

PAllF 

106.85366 

28.17669 

2013-02-08 

GAC109N31M1 

109.37725 

31.69805 

2012-01-29 

PB05M 

70.09451 

29.97223 

2013-02-09 

GAC072N28B1 

72.31558 

28.35083 

2012-01-31 

PB08M 

91.59353 

29.51241 

2013-02-09 

GAC094N27B1 

94.58662 

27.21015 

2012-01-31 

PB09B 

100.33707 

28.19664 

2013-02-10 

GAC110N25B1 

110.86880 

25.05053 

2012-02-14 

PA08F 

124.98256 

39.88470 

2013-02-14 

GAC053N32B1 

53.75501 

32.01677 

2012-02-16 

GAC.106N28.fi 

106.85366 

28.17669 

2013-02-14 

GAC114N32B1 

114.97536 

32.00973 

2012-03-13 

PA13B 

129.34213 

28.29440 

2013-02-15 

GAC076N30B1 

76.06069 

30.49459 

2012-03-14 

PB09B 

100.33707 

28.19664 

2013-02-18 

GAC053N32B1 

53.75501 

32.01677 

2012-10-03 

GAC080N33B101 

79.75012 

33.74839 

2013-03-04 

GAC.098N33.B1 

98.42774 

33.02404 

2012-10-05 

GAC054N25B1 

54.29599 

25.99110 

2013-03-04 

GAC.128N36.B1 

128.34054 

36.43643 

2012-10-05 

GAC067N27B1 

67.33264 

27.40422 

2013-03-04 

GAC.128N36.B2 

128.34054 

36.43643 

2012-10-05 

GAC089N24B1 

89.23383 

24.24969 

2013-03-05 

GAC078N26B1 

78.08962 

26.45461 

2012-10-06 

GAC085N33B1 

85.74273 

33.31455 

2013-03-05 

GAC108N24B1 

108.10991 

24.12859 

2012-10-06 

GAC085N33B2 

85.74273 

33.31455 

2013-03-06 

GAC089N28B1 

89.14070 

28.94227 

2012-10-07 

GAC081N30B1 

81.78447 

30.20860 

2013-03-06 

GAC089N28B2 

89.14070 

28.94227 

2012-10-07 

GAC081N30B2 

81.78447 

30.20860 

2013-03-07 

GAC082N33B1 

82.00390 

33.76370 

2012-10-12 

GAC056N24B1 

56.20090 

24.28947 

2013-04-03 

test_114N22_Bl 

114.41732 

22.34006 

2012-10-12 

GAC056N24F1 

56.20090 

24.28947 
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22 A. Rebass a-Mans erg as et al. 


Table A2. Names, coordinates, plate-spectrograph-fibre identifiers, XSTPS-GAC magnitudes and stellar parameters (effective temper¬ 
ature, surface gravity and mass) of the 92 DA LSS-GAC white dwarfs identified in this work. 


Jname 

RA 

[deg] 

DEC 

[deg] 

plateid 

spid fiberid 

9 

r 

i 

TeS 

[k] 

err 

log(s) 

[dex] 

err 

M 

[M©] 

err 

J062647.53+264552.1 

96.69805 

26.76446 

GAC.097N28.fi 

5 

220 

17.53 

17.64 

18.03 

27425 

562 

8.15 

0.10 

0.73 

0.06 

J070057.53+284310.1 

105.23969 

28.71946 

GAC.105N29.M1 

5 

129 

17.37 

16.98 

17.23 

13173 

462 

8.27 

0.16 

0.78 

0.11 

J074722.19+295015.3 

116.84244 

29.83757 

GAG.118N28.fi 

16 

27 

18.19 

17.83 

18.05 

15601 

1150 

7.78 

0.28 

0.50 

0.15 

J074850.71+301003.5 

117.21131 

30.16765 

GAG.118N28.fi 

16 

228 

18.27 

17.89 

18.11 

- 

- 

- 

- 

- 

- 

J075106.48+301727.0 

117.77698 

30.29082 

GAC.118N28.M1 

11 

169 

15.65 

15.92 

16.39 

32222 

401 

8.04 

0.10 

0.67 

0.06 

J074742.05+280945.6 

116.92519 

28.16266 

GAG.118N28.fi 

3 

60 

17.83 

17.43 

17.69 

13234 

286 

8.18 

0.09 

0.72 

0.06 

J075251.35+271513.9 

118.21396 

27.25385 

GAC.118N28.M1 

5 

171 

16.72 

16.73 

17.15 

23343 

737 

7.78 

0.11 

0.52 

0.06 

J034037.94+230304.2 

55.15809 

23.05117 

GAC056N24F1 

5 

224 

18.07 

18.35 

18.35 

- 

- 

- 

- 

- 

- 

J051002.11+231541.0 

77.50881 

23.26140 

GAC079N24B1 

10 

62 

14.70 

15.26 

15.46 

11134 

131 

7.89 

0.11 

0.56 

0.07 

J060027.42+234311.2 

90.11425 

23.71979 

GAG089N24M1 

8 

7 

16.98 

16.98 

16.98 

8205 

113 

7.99 

0.22 

0.59 

0.15 

J033633.53+234938.7 

54.13971 

23.82742 

GAC056N24F1 

10 

57 

17.63 

17.93 

17.96 

30419 

1165 

7.88 

0.30 

0.58 

0.16 

J073145.10+235352.9 

112.93792 

23.89802 

GAC110N25B1 

7 

144 

16.00 

16.22 

16.45 

11736 

372 

7.75 

0.27 

0.47 

0.16 

J033900.07+242510.5 

54.75028 

24.41957 

GAC056N24F1 

3 

53 

18.01 

18.49 

18.52 

17304 

953 

7.88 

0.22 

0.55 

0.12 

J055727.92+243558.8 

89.36633 

24.59967 

GAG086N24M1 

6 

56 

16.47 

16.97 

17.27 

29387 

567 

8.04 

0.13 

0.67 

0.08 

J080230.40+244922.6 

120.62665 

24.82295 

GAG120N25M1 

4 

119 

17.66 

17.75 

17.87 

9070 

222 

8.15 

0.34 

0.69 

0.21 

J054335.91+250410.8 

85.89961 

25.06966 

GAG084N26M1 

7 

122 

15.90 

16.44 

16.70 

68101 

4213 

7.37 

0.20 

0.51 

0.06 

J034422.26+251453.3 

56.09274 

25.24813 

GAC058N25B1 

14 

66 

15.39 

15.71 

15.85 

7931 

101 

7.74 

0.29 

0.45 

0.18 

J062159.50+252335.9 

95.49790 

25.39331 

GAG094N27M1 

7 

21 

17.56 

17.62 

17.70 

11728 

651 

8.25 

0.31 

0.76 

0.19 

J032817.13+252853.5 

52.07136 

25.48152 

GAG051N24M1 

9 

36 

16.57 

16.98 

17.15 

14132 

944 

8.00 

0.16 

0.61 

0.10 

J032854.06+252626.3 

52.22527 

25.44064 

GAG051N24M1 

9 

35 

18.06 

17.73 

17.54 

- 

- 

- 

- 

- 

- 

J054613.53+255031.7 

86.55636 

25.84214 

GAC.082N27.M1 

2 

249 

17.33 

17.62 

17.78 

- 

- 

- 

- 

- 

- 

J070950.16+255303.6 

107.45899 

25.88434 

GAC.106N28.M1 

1 

184 

17.10 

17.63 

17.94 

77295 

10207 

7.64 

0.42 

0.61 

0.16 

J053727.46+260611.3 

84.36441 

26.10313 

GAG084N26M1 

4 

135 

16.76 

17.17 

17.33 

15782 

488 

7.95 

0.11 

0.59 

0.06 

J064452.84+260947.7 

101.22017 

26.16326 

GAC.100N28.B1 

7 

168 

15.48 

15.98 

16.22 

16910 

671 

7.60 

0.16 

0.42 

0.07 

J071223.81+260933.4 

108.09919 

26.15928 

GAC.106N28.M1 

7 

226 

16.85 

17.24 

17.41 

14132 

645 

8.00 

0.15 

0.61 

0.09 

J055046.50+261220.3 

87.69377 

26.20563 

GAC.082N27.B1 

5 

18 

15.13 

15.64 

15.91 

21289 

814 

7.60 

0.14 

0.43 

0.06 

J063532.49+261958.6 

98.88537 

26.33295 

GAC.063022N281243.M1 

7 

221 

16.63 

17.17 

17.47 

34526 

928 

8.01 

0.18 

0.66 

0.11 

J063828.24+263359.6 

99.61767 

26.56656 

GAC.100N28.B1 

2 

52 

16.03 

15.60 

15.48 

- 

- 

- 

- 

- 

- 

J063856.01+263022.5 

99.73339 

26.50626 

GAC.100N28.B1 

2 

51 

14.85 

14.35 

14.22 

- 

- 

- 

- 

- 

- 

J025737.25+264047.9 

44.40520 

26.67997 

GAC.045N28.M1 

2 

143 

16.91 

17.00 

17.16 

- 

- 

- 

- 

- 

- 

J063810.91+264040.9 

99.54544 

26.67802 

GAC.100N28.B1 

2 

54 

16.32 

15.59 

15.28 

- 

- 

- 

- 

- 

- 

J063919.88+264102.6 

99.83283 

26.68405 

GAC.100N28.B1 

2 

56 

15.97 

15.58 

15.40 

- 

- 

- 

- 

- 

- 

J070716.11+263857.5 

106.81713 

26.64930 

GAC.106N28.M1 

1 

154 

17.63 

18.07 

18.20 

14899 

783 

7.29 

0.19 

0.31 

0.06 

J034605.50+264348.3 

56.52290 

26.73008 

GAG.034118N284209.fi 

7 

240 

17.54 

17.95 

18.26 

22550 

2127 

7.95 

0.35 

0.60 

0.21 

J070755.01+265103.0 

106.97921 

26.85082 

GAC.106N28.B1 

5 

12 

15.53 

15.86 

16.01 

15071 

952 

7.89 

0.20 

0.55 

0.12 

J070706.33+265756.7 

106.77636 

26.96574 

GAC.106N28.M1 

5 

138 

17.19 

17.63 

17.83 

16910 

564 

8.21 

0.14 

0.74 

0.09 

J064743.70+270906.2 

101.93207 

27.15172 

GAC.096N28.M1 

6 

25 

16.00 

16.59 

16.90 

41991 

1207 

8.21 

0.12 

0.79 

0.07 

J034409.58+271507.3 

56.03992 

27.25202 

PB03M 

5 

172 

16.02 

16.52 

16.75 

19868 

246 

7.78 

0.05 

0.51 

0.02 

J054447.48+272032.0 

86.19783 

27.34223 

GAC.082N27.M1 

3 

76 

17.08 

16.93 

16.90 

6459 

56 

7.75 

0.14 

0.45 

0.07 

J071520.95+273433.4 

108.83731 

27.57594 

GAC.106N28.M1 

6 

192 

17.93 

18.01 

18.12 

- 

- 

- 

- 

- 

- 

J061540.92+275202.0 

93.92048 

27.86723 

GAG094N27M1 

15 

229 

16.47 

16.91 

17.20 

46042 

1008 

8.82 

0.11 

1.13 

0.05 

J035123.16+280256.6 

57.84651 

28.04905 

GAG.034118N284209.fi 

6 

169 

17.75 

18.11 

18.51 

32222 

862 

8.90 

0.19 

1.16 

0.09 

J061000.28+281426.9 

92.50117 

28.24081 

test.094N28.Bl 

10 

131 

14.33 

14.74 

15.10 

18542 

62 

7.97 

0.02 

0.60 

0.01 

J042945.25+282224.2 

67.43856 

28.37338 

GAC.068N28.B1 

3 

248 

14.60 

15.21 

15.52 

39642 

699 

8.06 

0.09 

0.70 

0.05 

J074741.42+282104.9 

116.92258 

28.35137 

GAG.118N28.fi 

3 

54 

18.08 

18.40 

18.64 

13284 

2679 

8.10 

0.43 

0.67 

0.27 

J080748.88+282626.4 

121.95366 

28.44067 

GAC.080935N290534.M1 

5 

55 

17.79 

17.93 

18.02 

9356 

256 

7.97 

0.51 

0.58 

0.32 

J052147.24+283532.5 

80.44682 

28.59236 

GAC.080N28.M1 

4 

198 

17.50 

17.73 

17.86 

10780 

199 

8.29 

0.15 

0.79 

0.09 

J030406.41+285143.2 

46.02670 

28.86200 

GAC.045N28.M1 

4 

206 

16.82 

17.19 

17.39 

- 

- 

- 

- 

- 

- 

J030214.72+285707.4 

45.56134 

28.95206 

GAC.045N28.M1 

15 

183 

17.21 

17.60 

17.80 

18120 

1667 

7.51 

0.35 

0.39 

0.15 

J053931.86+285456.7 

84.88277 

28.91574 

GAC.079N29.M1 

6 

167 

17.39 

16.64 

16.17 

- 

- 

- 

- 

- 

- 

J042334.22+290205.7 

65.89260 

29.03492 

GAC.068N28.fi 

14 

44 

18.10 

18.32 

18.38 

24443 

2180 

8.06 

0.35 

0.67 

0.21 

J080800.00+290152.6 

122.00001 

29.03127 

GAC.080935N290534.M1 

4 

155 

16.64 

17.19 

17.48 

24163 

600 

7.78 

0.09 

0.52 

0.04 

J033134.40+291321.4 

52.89335 

29.22260 

GAC.055N28.M1 

14 

111 

17.14 

17.40 

17.55 

14393 

1093 

7.56 

0.26 

0.39 

0.12 

J040403.73+291703.8 

61.01555 

29.28440 

GAC.060N28.M1 

9 

165 

17.73 

17.80 

17.73 

- 

- 

- 

- 

- 

- 

J040842.17+292130.3 

62.17569 

29.35843 

GAC.063N29.M1 

10 

223 

17.06 

17.01 

16.97 

7659 

113 

7.95 

0.23 

0.56 

0.15 

J071004.83+292402.8 

107.52014 

29.40079 

GAC.107N27.B1 

16 

178 

15.29 

15.62 

15.87 

14089 

732 

7.94 

0.18 

0.57 

0.11 

J042435.27+293651.7 

66.14697 

29.61437 

GAC.068N28.M1 

14 

104 

17.29 

17.52 

17.69 

36572 

1719 

8.22 

0.30 

0.79 

0.18 

J054658.08+293633.4 

86.74200 

29.60928 

PA09M1 

14 

111 

17.23 

17.29 

17.16 

28064 

789 

7.84 

0.16 

0.56 

0.08 

J065740.34+300915.8 

104.41810 

30.15438 

GAC.105N29.M1 

3 

134 

17.14 

17.08 

17.06 

7839 

43 

7.90 

0.11 

0.54 

0.07 

J071037.96+301146.4 

107.65815 

30.19622 

GAC109N30M1 

10 

26 

17.23 

17.30 

17.38 

- 

- 

- 

- 

- 

- 

J030128.01+301536.6 

45.36671 

30.26016 

GAC.045N28.B1 

11 

19 

14.75 

15.18 

15.36 

- 

- 

- 

- 

- 

- 

J071011.40+303041.7 

107.54751 

30.51158 

GAC.106N28.M1 

11 

27 

16.46 

16.84 

17.10 

18756 

594 

8.33 

0.12 

0.82 

0.08 
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DA white dwarfs from the LSS-GAC 23 


Table A2. Continued. 


Jname 

RA 

[deg] 

DEC 

[deg] 

plateid 

spid fiberid 

9 

r 

i 

TeH 

[k] 

err 

log(9) 

[dex] 

err 

M 

[M©] 

err 

J051624.54+303910.3 

79.10227 

30.65286 

GAC080N32F1 

2 

118 

17.84 

18.14 

18.33 

21535 

2125 

8.46 

0.37 

0.91 

0.21 

J060113.75+304014.3 

90.30730 

30.67064 

PA09M1 

12 

66 

17.49 

17.66 

17.73 

10939 

330 

8.12 

0.31 

0.68 

0.20 

J052038.36+304822.6 

80.15984 

30.80629 

GAC.078N28.B1 

16 

139 

15.38 

15.68 

15.88 

15601 

395 

7.87 

0.09 

0.54 

0.05 

J033149.68+305944.9 

52.95702 

30.99581 

GAC.052N30.M1 

9 

224 

16.89 

17.41 

17.67 

19868 

306 

8.46 

0.06 

0.91 

0.04 

J080039.67+305537.3 

120.16529 

30.92703 

GAG121N33M1 

2 

118 

17.44 

17.84 

18.18 

25302 

1399 

7.58 

0.20 

0.44 

0.08 

J031430.60+310301.4 

48.62749 

31.05040 

GAG049N32M1 

5 

42 

17.20 

17.61 

17.82 

18120 

1669 

8.03 

0.38 

0.64 

0.23 

J072438.30+310729.0 

111.15957 

31.12471 

GAG109N30M1 

9 

238 

17.35 

17.51 

17.69 

10772 

381 

8.79 

0.31 

1.09 

0.11 

J075724.50+311310.3 

119.35207 

31.21953 

GAG121N33F1 

2 

147 

17.73 

18.19 

18.41 

20098 

741 

7.85 

0.14 

0.54 

0.07 

J063058.95+312344.1 

97.74564 

31.39559 

GAG096N32M1 

6 

239 

17.24 

17.53 

17.76 

15071 

803 

7.72 

0.19 

0.47 

0.10 

J055234.17+312401.6 

88.14237 

31.40045 

GAG090N33B1 

2 

89 

14.64 

15.00 

15.22 

14228 

176 

8.15 

0.03 

0.70 

0.02 

J082705.52+313008.2 

126.77302 

31.50228 

test_126N31_Bl 

9 

22 

15.36 

15.92 

16.32 

71310 

4547 

7.44 

0.22 

0.54 

0.06 

J071603.19+315711.0 

109.01331 

31.95305 

GAG109N30M1 

15 

11 

17.72 

17.91 

17.86 

- 

- 

- 

- 

- 

- 

J080602.45+315431.5 

121.51019 

31.90875 

GAG121N33M1 

8 

185 

17.25 

17.80 

18.11 

22811 

1761 

7.78 

0.28 

0.52 

0.15 

J062856.67+320303.9 

97.23611 

32.05109 

GAG098N33F1 

5 

193 

17.77 

18.31 

18.75 

36572 

1334 

7.71 

0.24 

0.53 

0.11 

J071451.13+320407.9 

108.71306 

32.06885 

GAG109N30M1 

15 

76 

15.89 

16.44 

16.79 

49906 

2195 

8.49 

0.18 

0.97 

0.09 

J053712.39+321502.9 

84.30162 

32.25081 

GAG085N31M1 

15 

243 

16.70 

16.96 

17.19 

16525 

275 

7.66 

0.06 

0.44 

0.03 

J042006.39+323305.0 

65.02661 

32.55138 

GAG065N31B1 

4 

216 

15.76 

16.20 

16.43 

17912 

544 

8.06 

0.11 

0.65 

0.07 

J063100.30+324453.8 

97.75127 

32.74827 

GAG100N32M1 

14 

38 

16.16 

16.76 

17.09 

64291 

7505 

7.70 

0.41 

0.60 

0.18 

J031336.68+325107.9 

48.40285 

32.85219 

GAG049N32M1 

15 

248 

17.29 

17.40 

17.52 

10642 

294 

8.43 

0.26 

0.87 

0.14 

J031448.30+324916.3 

48.70124 

32.82120 

GAG049N32M1 

15 

195 

16.70 

17.12 

17.32 

25595 

915 

8.24 

0.13 

0.78 

0.08 

J054709.65+324843.7 

86.79021 

32.81214 

GAG085N33M1 

8 

2 

16.72 

16.31 

16.13 

- 

- 

- 

- 

- 

- 

J061736.85+325732.0 

94.40353 

32.95888 

GAG096N32M1 

14 

173 

16.49 

16.97 

17.20 

35331 

576 

7.88 

0.12 

0.60 

0.06 

J080006.18+325738.1 

120.02574 

32.96058 

GAG121N33F1 

3 

179 

18.15 

18.45 

18.63 

14393 

1410 

8.44 

0.20 

0.89 

0.13 

J030945.71+330025.0 

47.44047 

33.00694 

GAG049N32M1 

14 

79 

16.78 

17.09 

17.30 

14899 

1087 

8.62 

0.18 

1.00 

0.10 

J075321.55+334308.1 

118.33978 

33.71892 

GAG121N33F1 

14 

103 

17.98 

18.41 

18.68 

36996 

2504 

7.07 

0.41 

0.34 

0.11 

J055407.71+340348.0 

88.53214 

34.06334 

GAG087N32M1 

11 

71 

17.59 

17.62 

17.76 

10409 

188 

7.19 

0.40 

0.27 

0.17 

J055019.78+350006.5 

87.58242 

35.00180 

GAG085N33M1 

12 

184 

16.31 

16.64 

16.91 

18542 

404 

7.70 

0.08 

0.47 

0.04 

J051940.72+355137.1 

79.91965 

35.86031 

GAG080N33M1 

11 

178 

17.27 

17.45 

17.68 

10656 

174 

8.08 

0.21 

0.65 

0.14 

J083400.13+365921.7 

128.50056 

36.98935 

GAG.128N36.B2 

4 

206 

16.05 

16.09 

16.14 

9065 

67 

7.93 

0.12 

0.56 

0.08 

J084218.71+374859.4 

130.57797 

37.81649 

GAG.128N36.B2 

12 

116 

16.05 

16.46 

16.76 

22037 

1183 

8.42 

0.20 

0.88 

0.12 
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